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Abstract. We considerably improve upon the recent result of [34] on the mixing 
time of Glauber dynamics for the 2d Ising model in a box of side L at low temperature 
and with random boundary conditions whose distribution P stochastically dominates 
the extremal plus phase. An important special case is when P is concentrated on the 
homogeneous all-plus configuration, where the mixing time Tmix is conjectured to be 
polynomial in L. In 34 it was shown that for a large enough inverse-temperature /3 
and any £ > there exists c — c{l3, s) such that Huil^oo P (Tmix > exp(cL^)) — 0. In 
particular, for the all-plus boundary conditions and /3 large enough Tmix < exp(cL^). 

Here we show that the same conclusions hold for all /3 larger than the critical 
value /3c and with exp(cL^) replaced by L^^^^^ (^.g, quasi-polynomial mixing). The 
key point is a modification of the inductive scheme of 34 together with refined 
equilibrium estimates that hold up to criticality, obtained via duality and random- 
line representation tools for the Ising model. In particular, we establish new precise 
bounds on the law of Peierls contours which quantitatively sharpen the Brownian 
bridge picture established e.g. in |17lp!9|[20] . 

1. Introduction 

The Ising model on lattices at and near criticality has been the focus of numerous 
research papers since its introduction in 1925, establishing it as one of the most studied 
models in mathematical physics. In two dimensions the model was exactly solved by 
Onsager |37| in 1944, determining its critical inverse-temperature 13c — \ log(l + a/2) in 
the absence of an external magnetic field. While the classical study of the Ising model 
concentrated on its static properties, over the last three decades significant efforts were 
dedicated to the analysis of stochastic dynamical systems that both model its evolution 
and provide efficient methods of sampling it. Of particular interest is the interplay 
between the behaviors of the static and dynamical models as they both undergo a 
phase transition at the critical /3c. 

The Glauber dynamics for the Ising model (also known as the stochastic Ising model) ^ 



introduced by Glauber 16 in 1963, is considered to be the most natural sampling 
method for it, with notable examples including heat-hath and Metropolis. It is known 
that on a box of side-length L in I? with free boundary conditions (b.c), alongside 
the phase transition in the range of spin-spin correlations in the static Ising model 
around /3c, the corresponding Glauber dynamics exhibits a critical slowdown: Its mixing 



time (formally defined in ^1.1) transitions from being logarithmic in L in the high 
temperature regime /3 < /3c to being exponentially large in L in the low temperature 
regime ^5 > /3c, en route following a power law at the critical /3c. 
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One of the most fundamental open problems in the study of the stochastic Ising 
model is understanding the system's behavior in the so-called phase-coexistence re- 
gion under homogenous boundary conditions, e.g. all-plus boundary. In the presence 
of these b.c. the (— ) phase becomes unstable and as such the reduced bottleneck be- 
tween the two phases drastically accelerates the rate of convergence of the dynamics 
to equilibrium. Indeed, in this case the Glauber dynamics is known to mix in time 
that is sub-exponential in the surface area of the box, contrary to its low-temperature 
behavior with free boundary. The central and longstanding conjecture addressing this 
phenomenon states that the mixing time of Glauber dynamics for the Ising model on 
a box of side-length L with all-plus boundary conditions is at most polynomial in L at 
any temperature. 

So far this has been confirmed on the 2d lattice throughout the one-phase region 
/3 < /3c (see [30j|3l]) and very recently at the critical (3 — 13c (see [27]). Despite 
intensive efforts over the last two decades, establishing a power-law behavior for the 
mixing of Glauber dynamics at the phase-coexistence region /3 > /3c under the all-plus 
b.c. remains an enticing open problem. 

In the precise order of mixing in this regime on a 2d square lattice of side-length L 



36 



44 



was conjectured to be L? in accordance with Lifshitz's law (see [25] and also [9 
The heuristic behind this prediction argues that when a droplet of the (— ) phase is 
surrounded by the (+) phase at low temperature it proceeds to shrink according to the 
mean-curvature of the interface between them. Unfortunately, rigorous analysis is still 
quite far from establishing the expected Lifshitz behavior of 0{L'^) mixing. 

Until recently the best upper bound on the mixing at the phase-coexistence region 
under the all-plus boundary was exp(L-'^/^+^(-'^)) due to |29| and valid for large enough /3. 
This bound from 1994 was substantially improved in a recent breakthrough paper |34 , 
where it was shown (as a special case of a result on a wider class of b.c.) that for 
a sufficiently large /3 and any e > the mixing time is exp(0(L^)). The approach 
of 1 34 1 hinged on a novel inductive scheme on boxes with random boundary conditions, 
combined with a careful use of the so-called Peres- Winkler censoring inequality] these 
ideas form the foundation of the present paper. Note that the requirement of large /3 
in 29,34) was essential in order to make use of results of [IT] on the Wulff construction. 



available only at low enough temperature by cluster expansion methods. For smaller 
values of /3 > /3c the best known estimates on the mixing time are due to |8| and of the 
weaker form exp(o(L)). 

In this work we improve these estimates into an upper bound of the form L^^^^^^^ 
on the mixing-time (i.e. quasi-polynomial in the side-length L) valid for any /3 > Pc- 
The key to our analysis is a modification of the recursive framework introduced in [34j 
combined with refined equilibrium estimates that hold up to criticality. To establish 
these, in lieu of relying on cluster-expansions, we utilize duality and the random-line 
representation machinery for the high temperature Ising model developed in [40{[4T] . 

A key new element of our proof concerns fine estimates on the fluctuations of cluster 
boundaries. Whenever the boundary is pinned at two vertices u = (0, 0) and v = (£, 0), 
the contour of the cluster is known to converge to the Brownian bridge (cf. [To|[T9l[20] ). 
This does not, however, allow us to estimate the probability of events when these 
converge to in the limit. In particular, we are interested in: (i) a Gaussian boun d fo r 
the probability that the contour would reach height xVI (established in Theorem 5.3); 
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(ii) the probabihty that the contour remains in the upper half-plane, an event that 
would have probability l/l were the contour to behave like a Id random walk of length 
I conditioned to return to 0. In ^ (see Theorem 5.1 ) we prove that up to multiplicative 
constants this indeed holds for a given contour. 

These then provide important tools in estimating the probability of various other 
events characterizing the Ising interfaces at equilibrium. 

1.1. Glauber dynamics for the Ising model. Let A be a generic finite subset of I?. 
Write X ^ y for the nearest-neighbor relation in [i.e. x ^ y if max^^i^2 \xi — yi\ — 1) 
and define 5A, the boundary of A, to be the nearest-neig hbors of A in \ A: 

5A := {x G \ A : x ^ y for some y G A} . 

The classical Ising model on A with no external magnetic field is a spin-system whose set 
of possible configurations is I^a = { — 1, +1}^. Each configuration a G f^A corresponds 
to an assignment of plus/minus spins to the sites in A and has a statistical weight 
determined by the Hamiltonian 

x.yeA xeA.yedA 
x^y x^y 

where r G fi^A forms the boundary conditions (b.c.) of the system. The Gibbs measure 
associated to the spin-system with boundary conditions r is 



^X(^) = ^e-/^""^^-) (aGf^A), (1.1) 



^^,A 

where /3 is the inverse of the temperature {i.e. /3 = ^) and the partition-function 
is a normalizing constant. When the boundary conditions are uniformly equal to +1 
(resp. —1) we will denote the Gibbs measure by tt^ (resp. tt^). Throughout the paper 
we will omit the superscript r and the subscript A from the notation of the Gibbs 
measure tt when these are made clear from the context. 

The Gibbs measure enjoys a useful monotonicity property that will play a key role 
in our analysis. Consider the usual partial order on Q\ whereby cr < r] if ax < r]x for 
all X G A. A function / : Q\ i-^ M is monotone increasing (decreasing) if cr < 77 implies 
f{o-) < f{r]) {f{cr) > f{r])). An event is increasing (decreasing) if its characteristic 
function is increasing (decreasing). Given two probability measures fi, v on Vl^ we say 
that [1 is stochastically dominated by z/, denoted by ^ z/, if ii{f) < Ty{f) for all 
increasing functions / (here and in what follows fi{f) stands for J fdji). According to 
these notations the well-known FKG inequalities [14J state that 

• If T < then ti\ ^ tt]^^ . 

• If / and g are increasing then Ti\{fg) > Ti\{f)Ti\{g) . 

The phase transition regime in the 2d Ising model occurs at low temperature and it 
is characterized by spontaneous magnetization in the thermodynamic limit. There is a 
critical value /3c such that for all /3 > /3c, 

1™, ^A (^0) = - lim 7r;^(cro) ^m^>{). (1.2) 
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Furthermore, in the thermodynamic hmit the measures tt^ and tt^ converge (weakly) 
to two distinct Gibbs measures tt^ and tt^ which are measures on the space , each 
representing a pure state. We wih focus on this phase-coexistence region j3 > j3c. 

The Glauber dynamics for the Ising model is a family of continuous-time Markov 
chains on the state space I^a, reversible with respect to the Gibbs distribution tt^. An 
important and natural example of this stochastic dynamics is the heat-bath dynamics, 
which we will now describe, postponing the formulation of the general Glauber dynam- 



ics to ^2.1 Note that our results apply to all of these chains (e.g.. Metropolis etc.) by 
standard arguments for comparing their mixing times (see e.g. p8]). 

The heat-hath dynamics for the Ising model Vt^ is defined as follows. With a rate 
one independent Poisson process for each vertex x, the spin cjx is refreshed by sampling 
a new value from the set { — 1, +1} according to the conditional Gibbs measure 

•= ^A (• \(^y^y^x) . 

It is easy to verify that the heat-bath chain is indeed reversible with respect to ti\ and 
is characterized by the generator 



xGA 



where TT^^xif) average of / with respect to the conditional Gibbs measure tt^^^ 

acting only on the variable ax- The Dirichlet form associated to C\ takes the form 

W,/) = E^A(Var;.(/)) 

xeA 

where VarJ .p(/) denotes the variance with respect to vrj ^,. It is possible to extend the 
above definition of the generator C\ directly to the whole lattice and get a well 
defined Markov process on fi^2 (see e.g. |26|). The latter will be referred to as the 
infinite volume Glauber dynamics, with generator denoted by C^2. 

We will denote by /i^ the distribution of the chain at time t when the starting 
configuration is identically equal to a. For instance, for any / and a the expectation of 
/ w.r.t. /if is given by (T]^(t))f(a) where TJ^ is the Markov semigroup generated by C\. 
The notation /x^ will stand for the corresponding quantity for an initial configuration 
of either all-plus or all-minus. 

A key quantity that measures the rate of convergence of Glauber dynamics to sta- 
tionarity is the gap in the spectrum of its generator, denoted by gap]^. The Dirichlet 
form associated with C\ produces the following characterization for the spectral-gap: 



gaPA = 



£lifJ) 



f VarX(/) ' 

where the infimum is over all nonconstant / G L'^{7tJJ. Another useful measure for the 
speed of relaxation to equilibrium is the total- variation mixing time which is defined as 
follows. Recall that the total- variation distance between two measures 0, '0 on a finite 
probability space is defined as 

U - sup \ct>{A) - i,iA)\ ^IYI - • 
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For any e G (0, 1), the 6-mixing-time of the Glauber dynamics is given by 
Tmix(£:) := infjt > : sup - ttXH < e\. 



When e = l/(2e) we wih simply write Tmix- This particular definition yields the 
following well-known inequalities (see e.g. |24l|42]): 



sup ll/i^ — ttXII < exp(— [^/TmixJ) for alH > , 

1 1 2e 

< Tmix < log where TTmin min 7rX(cr) . 

gap gap TTmin ^^^^A 

The last inequality shows that in our setting gap~^ and Tmix are always within a factor 
of 0(|A|) from one another (to see this, observe that Ti\{a) / Ti\{a') < exp(0(|A|)) for 
any a, G ft\ by Eq. ( |1.1| ) whereas |riA| = 2l^l). One could restate our results as well 
as the analogous conjecture on the polynomial mixing time under all-plus b.c. in terms 
of gap~^ (expected to have order L, the side-length of A, for any /3 > /3c; see [5||7J). 

1.2. Main results. We are now in a position to formalize the main contribution of 
this paper. The following theorem is the counterpart of the main result obtained by 
two of the authors in |34|. Here we feature an improved estimate that in addition holds 
not only for large enough /3 but throughout the phase-coexistence region. 

Theorem 1. For any /3 > /3c there exists some c(/3) > so that the following holds for 
the Glauber dynamics for the Ising model on the square Al at inverse-temperature /3. 
// L is of the form L — 2^ — 1 for some integer n then: 

(1) If the boundary conditions r are sampled from a law P that either stochastically 
dominates the pure phase tt^ or is stochastically dominated by tt^ then 

m^t - ^^11 ^ ^/^ f^^ = L''^''^^ . (1.3) 

In particular, 

P{Tu^x>L'^'°^'')<c/L. (1.4) 



(2) The estimates (1.3);( |1.4[ ) also hold ifF is stochastically dominated by tt^ on 
one side of and stochastically dominates tt^ on the union of the other three 
sides. A similar statement holds if the roles of (+) and (— ) are reversed. 

The most natural consequence of the above result is obtained when P concentrates 
on homogenous boundary conditions, where the best previous bounds were exp(0(L^)) 
for any e > and /3 large enough (|34|) along with exp(o(L)) for all other ^ > (|8|). 

Corollary 2. For any /3 > /3c there exists some c(/3) > 50 that the mixing time of 
Glauber dynamics for the Ising model on the square Al with b.c. r = +1 satisfies 

rM,x<i:^'°^''. (1.5) 

The same bound holds if the boundary conditions are (+) on three sides and (— ) on the 
remaining one, and similarly if (+) is replaced by (— ). 
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We believe that improving the above L^^^^^ bound into the conjectured polynomial 
one would require substantial new ideas. Indeed, in the present recursive framework 
in which the final scale of the system is reached via a doubling sequence, at each step 
the mixing-time estimate worsens by a power of L (hence the quasi-polynomial bound). 
For a polynomial upper bound one could not afford to lose more than a constant factor 
on average along these steps. 

One may also apply Theorem [l] to deduce the mixing behavior of the 2d Ising model 
under Bernoulli boundary conditions, as illustrated by the next corollary. Here and in 
what follows we say that an event holds with high probability (w.h.p.) to denote that 
its probability tends to 1 as the size of the system tends to oc. 

Corollary 3. Let /3 > /3c and consider Glauber dynamics for the Ising model on the 
square Al with b.c. r comprised of i.i.d. Bernoulli variables, P(t(x) = +1) = for 
some p+ > ^(1 + tanh(4/3)). Then w.h.p. T^ix < L^^^^^ for some c{P) > 0. 

To obtain the above corollary observe that the Bernoulli boundary conditions with 
the above specified p^ clearly stochastically dominate the marginal of tt^ on dA. 

The mixing time of Glauber dynamics for Ising on a finite box under all-plus b.c. 
is closely related to the asymptotic decay of the time auto-correlation function in the 
infinite-volume dynamics on started at the plus phase. Here it was conjectured 
that the decay should follow a stretched exponential of the form exp(— c\/t). As 



m 



12j 

a by-product of Corollary [2] (and standard monotonicity arguments) we obtain a new 
bound on this quantity, improving on the previous estimate due to |34| of exp( — (log t)*^) 
with arbitrarily large c which was applicable for large enough /3. 

Corollary 4. Let (3 > let f{a) = ctq and define p{t) = Var^ (e^^^^ j) iq ^/^g time 
autocorrelation of the spin at the origin started from the plus phase tt^ (the variance 
is w.r.t. the plus phase vr^^. Then there exists some c(/3) > such that for any t, 

p{t) <exp(-e^^^). (1.6) 

1.3. Related work. Over the last two decades considerable effort was devoted to the 
formidable problem of establishing polynomial mixing for the stochastic Ising model on 
a finite lattice with all-plus b.c. Following is a partial account of related results. 

Analogous to its conjectured behavior on Z^, the mixing of Glauber dynamics for the 
Ising model on the lattice in any fixed dimension d is believed to be polynomial in 
the side-length of the box at any temperature in the presence of an all-plus boundary. 
Unfortunately, the state-of-the-art rigorous analysis of the problem in three dimensions 
and higher is far more limited. Faced with the polynomial lower bounds of |5 , the 
best known upper bound for dimension d > 3 is exp(L^~^^^(-'^^) for large enough j3 (as 
usual L being the side-length) due to |45j. Compare this with the case of no {i.e. free) 
boundary conditions case where it was shown in [46j that gap~^ (and thus also Tmix) 
is at least cexp(c'^L^~-'^) for some c = c{l3) > and an absolute constant > 0. 



In two dimensions, ever since the work of Martinelli 29 in 1994 (an upper bound 
of exp(LV2+o(i)) at low enough temperatures) and until quite recently no real progress 
has been made on the original problem. Nevertheless, various variants of this problem 
became fairly well understood. For instance, nearly homogenous boundary conditions 
were studied in |2,3 . Analogues of the problem on non-amenable geometries (in terms 
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of a suitable parameter measuring the growth of bahs to replace the side-length) were 
estabhshed, pioneered by the work of ^ on trees and followed by results of |4 on a class 
of hyperbolic graphs of large degrees. The Solid-On-Solid model (SOS), proposed as an 
idealization of the behavior of Ising contours at low temperatures, was studied in |32 
where the authors obtained several insights into the evolution of the contours. Finally, 



the conjectured Lifshitz behavior of O(L^) was confirmed at zero temperature |Q^,13 
with the recent work |7| providing sharp bounds also for near-zero temperatures (namely 
when /3 > clogL for a suitably large c > 0) in both dimensions two and three. 

As mentioned above, the exp(L^/^) barrier was finally broken in the recent paper |34 , 
replacing it by exp(cL^) for an arbitrarily small e > and sufficiently large /3 (where 
the constant c = c(/3,6:) diverges to +oc as 5 ^ 0). At the heart of the proof of the 
main result of that paper (|34, Theorem 1.6]) was an inductive procedure which will 
serve as our main benchmark here. We will shortly review that argument in ^in order 
to motivate and better understand the new steps gained in the present work. 

Finally, there is an extensive literature on the phase-separation lines in the 2d Ising 
model, going back to In ^ we will review the tools we will need from the 

random-line representation framework of 40, 4l]. For further information see e.g. 39 
and the references therein. 

2. PRELIMINARIES 

2.1. General Glauber dynamics. The class of Glauber dynamics for the Ising model 
on a finite box A C consists of the continuous-time Markov chains on the state space 
f^A that are given by the generator 

(A/)(c7) = Yl ^) (/('^') - /(^)) ' (2-1) 

where is the configuration a with the spin at x flipped and the transition rates 
c(x, a) should satisfy the following conditions: 

(1) Finite range interactions: For some fixed R > and any x G A, if a, G f^A agree 
on the ball of diameter R about x then c(x, a) = c(x, a^). 

(2) Detailed balance: For all a G I^a and x G A, 

c(x, a) 



c(x, 



= exp(-/3Va.ifX(<7)), 



where VxHl{a) = Hl{a^) - HJ^{a) = 2ax J2y<^A cry + J2y€dA Ty 

(3) Positivity and houndedness: The rates c(x,a) are uniformly bounded from below 
and above by some fixed Ci, C2 > 0. 

(4) Translation invariance: If a = cr'(- + ^), where ^ G A and addition is according to 
the lattice metric, then c(x, a) = c(x + ^, a') for all x G A. 

The Glauber dynamics generator with such rates defines a unique Markov process, 
reversible with respect to the Gibbs measure fi^. The two most notable examples for 
the choice of transition rates are 

(i) Metropolis: c{x,a) = ex.p(^—^VxH]^{cr)^ A 1. 

(ii) Heat-bath: c(x,cr)= 1 + exp(/3Va:i?A(^))' 
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See e.g. |28| for standard comparisons between these chains, in particular implying that 
their individual mixing times are within a factor of at most 0(|A|) from one another 
(hence our results apply to every one of these chains). 

2.2. Surface tension. Denote by T/3(0) the surface tension that corresponds to the 
angle defined as follows. Associate with each angle 9 G [— 7r/4, 7r/4] the unit vector 
= (cos6>,sin6>) G and the following b.c. r]e for Al = {-[L/2J, . . . , [L/2J}2: 

. . / +1 if {v,ne) > 0, 

Let be the partition-function of the corresponding Ising model and, as usual, let 

del 

direction orthogonal to no is the limit 



^^Al denote the partition-function under the all-plus b.c. The surface tension in the 



= - lini log , 

which gives rise to an even analytic function with period 7r/2 on R (a closed formula 
appears e.g. in (io', Section 5]). One can then extend the definition of to M? by 
homogeneity, setting T/3(x) = ri3{0x)\x\, where |x| denotes the Euclidean norm of x and 
0x is the angle it forms with (1,0). For all j3 > j3c this qualifies as a norm on 

The surface tension measures the effect of the interface induced by the boundary 
conditions rfe on the free-energy and thus plays an important role in the geometry of 
the low temperature Ising model. For instance, it was shown in [43] that the large 
deviations of the magnetization in a square are governed by t^(0) (also see [2T||22]). 



One of the useful properties of the surface tension is the sharp triangle inequality 
(see for instance ^] Proposition 2.1]): For any /3 > /3c there exists a strictly positive 
constant t<i^ such that for any x, y G we have 

r(3{x) + Ti3{y) - T(3{x + y)> 1^13 {\x\ + \y\ - |x + y\) , (2.2) 

A thorough account of additional properties of the surface tension may be found e.g. 
in fn] and [39]. 

2.3. Duality. Let Z^* := + (|, |) denote the dual lattice to Z^. The collection of 
edges of Z^ and of Z^* will be denoted by E^2 and E^2* respectively. It is useful to 
identify an edge e = {x,y) G E^2 with the closed unit segment in M? whose endpoints 
are and similarly do so for edges in E^2*- To each edge e = (x^y) G E^2 there 

corresponds a unique dual edge e* G £"^2* defined by the condition e H e* 7^ 0. 

Given a finite box A C Z^ of the form A = {x = (xi,X2) G Z^ : \xi\ < Li, \x2\ < i^2}, 
the dual box A* C Z^* is A* = {x = (xi,X2) G Z^* : \xi\ < Li + ^, \x2\ < ^2 + 
The set of dual edges of A*, denoted by E"^*, is the set of dual edges for which both 
endpoints lie in A*. Notice that for each edge e = (x, y) G E^2 such that {x, H A 7^ 0, 
the corresponding dual edge e* necessarily belongs to These definitions readily 

generalize to an arbitrary finite ^ C Z^, in which case C Z^* consists of all dual 
sites whose L^-distance from Q equals 1. 

For any (3 > /3c we associate the dual inverse-temperature /3* via the duality relation 
tanh(/3*) = e~^^. Notice that for any /3 > /3c the dual inverse temperature /3* lies below 
/3c which is the unique fixed point of the map /3 ^5*. We will often refer to the Gibbs 
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Figure 1 . SE and S W splitting-rules for forming the contours 



measure on a subset of the dual lattice A* at the inverse-temperature /3* under free 
boundary, denoting it by tt^* . The following well-known fact addresses the exponential 
decay of the two-point correlation function for the free Ising Gibbs measure above the 
critical temperature. 



Lemma 2.1 (e.g. [35, p309 Eq. (4.39)], together with the GKS inequalities [T8|[23)). 
Let K dl? and /3 > /3c. There exists some C/3 > such that for any x, ^ G A*, 

^A* {^x (^y) < ( /— ^ A 1 ) exp ( - ri3{x - y)) . 

A matching exponent for the spin-spin correlation was established by [Tt] for two 
opposite points in the (dual) infinite strip. Let 5 = {1, . . . , ^} x Z for some integer £ 
and fix /3 > /3c. In the dual 5* we let x = (i, i) and y = (£+ i, i) and consider the free 



Gibbs measure at inverse-temperature /3*. It was shown in [T7| formula (2.22)] that in 
this setting there exists some > such that 

C/5 + 0(1) 



TT^-* (CT, 



■exp(-r^(0)£). 



(2.3) 



OC. 



where the o(l)-term tends to as ^ 

2.4. Contours. Let Q = {V,E) be a finite subgraph of (Z^*,£;*2*). The boundary of 
a subset of dual edges B d denoted by 5B^ is the set of vertices of V with an odd 
number of adjacent edges of 5. If 5B = we say that B is closed^ otherwise it is open. 

A chain of sites of length k from x to ^ in ^ has the standard definition of a sequence 
of sites X — uo^ui^ . . . — y such that ui ^ Q and \ui — Ui-i\ — 1 for all i. A ^- chain 
from X to y is similarly defined with the exception that the distance requirement is 
relaxed into 1 <\ui — Ui-i\ < \/2 for all z. A path from x to y in S is a chain of sites 
consisting of edges of that is {ui-\^ Ui) ^ B for all i. We say that a path is closed if 
its endpoint and starting point coincide, otherwise we say that it is open. 

A set of dual edges B C E can be uniquely partitioned into a finite number of 
edge-disjoint simple lines in H called contours. This is achieved by repeating the 
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Figure 2. Contour decomposition of the edge-set induced by an Ising 
configuration on a box with mixed boundary conditions according to 
both SE and SW sphtting-rules. 

following procedure referred to as the South- East (SE) splitting-rule: When four bonds 
meet at a vertex we separate them along the SE-oriented diagonal going through the 
intersection. Alternatively, one may globally apply the SW splitting-rule, analogously 
defined with the South- West orientation replacing the South-East one (see Figure [T]). 

Contours can be either open or closed (with the same distinction as in paths). The 
length of a contour 7, denoted by I7I, is the number of edges in 7, and the length of 
a collection of contours 7 = {71, . . . ,7^^} will simply be the sum of all the individual 
lengths. Given a finite family of contours 7 = {71, . . . , 7^^^^} we say that it is compatible 
if it is the contour decomposition of its collection of dual edges Uiji. We further say 
that 7 is E -compatible (or Q -compatible) to emphasize that in addition all the edges of 
Uiji belong to E, the edge-set of Q. 

Given boundary conditions r G { — 1,1}^ and a box A, each spin-configuration a 
compatible with r outside A (i.e. ax = Tx for any x ^ A) can be uniquely specified by 
giving all the edges e = (x, y) G £^^2 such that ax ^ ay and {x, ^} H A ^ (that is, all 
edges whose endpoint sites disagree). Equivalently, one can specify the corresponding 
dual edges of A*. By applying the above contour decomposition we see that each 
configuration a compatible with r is uniquely characterized by its collection of closed 
and open contours (see Figure [2] for an illustration). The open contours obtained in 
this manner are called the phase- separation lines. 

It is clear that the boundary SX of the open contours belongs to 5A* and must 
coincide with a certain set V{r) uniquely specified by the boundary conditions r 
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(i.e. independent of the values a gives to the spins of A). Notice that the cardinal- 
ity of V{r), if different from zero, must be even. 

A family of closed and open simple lines 7 is called r- compatible if there exists a 
configuration a compatible with r in A from which 7 is obtained in the above procedure. 
One can easily verify that when A is a box the set of r-compatible contours coincides 
with the set of E"^* -compatible contours whose boundary is equal to V{r). 

2.5. Random-line representation. For a finite subgraph Q = {V,E) of (Z^*,£^^2*) 

and an E'-compatible family of contours 9, two different partition functions Z{G) and 
Z{Q I 0) will turn out to be useful for a given /3 > 0: 

z{g)= J2 (2.4) 

7: ^7=0 
7 is £^-compatible 

z{g\e)= Yl (2.5) 

7: ^7=0 
jUO is £^-compatible 

Using Z(Q) and Z(Q \ 9) we define the weight (not necessarily a probability distribution) 
corresponding to the family of contours 9, denoted by qg(9), to be 

/m f ^^e-2/5|^l if is £;-compatible, 
qg{0) = < /(^) - . (2.6) 

I otherwise. 

The key reason for the above formula is the following random-line representation for 
even-point correlation functions: Consider the Ising model on Q at inverse temperature 
/3* and free boundary conditions. Let ttq be the associated Gibbs measure and let 



A cV have even cardinality. Then the following holds (see 40, Lemma 6.9]): 



E ^e(A) = 7ra(n^x). (2.7) 

A: SX=A xeA 



Remark. If the cardinality of A is odd then the r.h.s. of (2.7) is zero by symmetry 
and the l.h.s. is zero due to the definition of qg{X). 

Back to the low temperature Ising model in a box A with boundary condition r, let 
A be a collection of r-compatible open contours. Then, by construction, 

^rf^. W^A ^A* (A) ^A* (A) .^^x 

7Tj^[a : A(cr) = Aj = = —7- . ^ , (2.8) 

Z^X':SX'=V{t) ^A* (A j TT^* ( I IxEVir) ^ x) 

where with a shght abuse of notation we have identified A* with the graph Q = (A*, £^^*) 



and in the last equality we used (2.7). The above formula will be the starting point of 



the proof of the new equilibrium estimates. Propositions 4.4 and 4.5, 

We conclude this section with some of the main properties of the weights qg{X). For 
further information see |40||4T]. 



Lemma 2.2 ( 40, Lemma 6.3]). Let Q = {V,E) he a finite subgraph 0/ (Z^*, £;*2*) and 
let 9^ he a family of E- compatible contours (open and closed). If Q' is a subgraph of Q 
then qgi (9) > qg (9) . 
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Remark. Lemma 2.2 enables one to extend the definition of the weights qg{X) for 



finite contours A in a infinite graph G by taking a limit for n ^ oc of qg^{X), where 



Gn is the intersection of G with a box of size n. By Lemma 2.2 the sequence qg^{X) is 
monotone non-increasing and non-negative hence its limit indeed exists. 

Let G = {V^E) be a subgraph of (Z^*,£^*2*). The edge-boundary of an edge e ^ 
denoted by A(e), is comprised of the edge e itself together with any edge e' ^ E that 
is incident to it and would belong to the same contour in the contour decomposition of 
E via the agreed splitting-rule. For instance, with the SE splitting-rule the horizontal 
edge e = [{x^y)^ {x + 1,^)] in the dual lattice Z^* would have an edge-boundary of 
A(e) = {e, [(x, y), (x, y + 1)], [{x + 1, ?/), (x + 1, — 1)]}. Given a subset of edges B C E 
we define its edge-boundary as A(S) = UeGBA(e). This definition implies that two 
contours A and 7, where A is closed and 7 is either open or closed, are ^-compatible if 
and only if the edge-set of A does not intersect A(7) (see the related [40[ Lemma 6.1]). 
The following lemma is a special case of Lemma 6.4]): 



Lemma 2.3 (|40| Eq. (6.17)]). Let G = {V,E) be a subgraph of (Z^*,^;*^,) and let 
9 and A denote two G -compatible families of contours with corresponding edge-sets Eq 
and Ex respectively. If XU 6 is G -compatible (or equivalently if A{X) D Eq^ then 

qg(9U X) = qg^(9)qg(X) , 

where G\ is the subgraph of G given by the edge-set E\A(X). 

We will frequently need estimates on the weight of a contour constrained to go 
through certain dual sites; to this end, the following definition will be useful. Let 
G = (V^ E) and let Ai, A2 be two open contours such that 5Xi — {x, y} and 8X2 — {u^ v}. 
We say that Ai,A2 are disjoint if either they are ^-compatible or their edge-sets are 
disjoint and the contour decomposition of the union of their edges is a single contour 
A. Observe that in the latter case necessarily {x, y} H {u^ v} ^ 0. For a pair of disjoint 
open contours Ai, A2 we write Ai U A2 to denote either the collection (Ai, A2) in the 
former case or the single contour A in the latter. 

Lemma 2.4 ([40( Lemma 6.5]). Let G = (^, E) be a graph in the dual lattice Z^*. For 
any x, y^u^v ^ V, 

A=AiUA2 Ai A2 

SXi={x,y} ,SX2={u,v} SXi={x,y} SX2={u,v} 

In particular. 

Corollary 2.5 (fir, Eq. (5.29)]). Let G = {V,E) be a graph in the dual lattice Z^*. 
For any /3 > /3c and any u,v,z G V, 

^qW^ [ ) ( ^^(^) ) = 4 (^^ i^v (yz) • 

X:dX={u,v} ^X:dX={u,z} ' ^X:6X={z,v} ' 

zex 



Together with Lemma 2.1 the above lemma immediately implies an upper bound 
on the weights in mention in terms of the surface tensions ri^(u — v) and ti^(v — z)). 
The next lemma provides an analogous bound for the weights of closed contours going 
through a set of prescribed sites. 
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Lemma 2.6 (|41, Lemma 5.5 part (ii)]). Let Q = {V,E) he a graph in Z^*. Let 
xi, . . . , x/e G y and identify xq = Xk- Then 



qg{X) <expf -^T;3(xi-Xi_i) j 

:SX=^ ^ i=l ' 



xi,...,a:fcGA 



3. Inductive framework for rectangles with "plus" boundaries 

In this section we outline the recursive scheme developed in |34 which, as mentioned 
in ^ established a significantly improved upper bound of exp(cL^) for the mixing time 
on a box of side-length L with "plus" b.c. at sufficiently low temperatures. 

Given 6 > (to be thought of as very small) and L G N let 

Rl-{x^ (ij) G : 1 < i < L, 1 < j < [Li+n}. 

Similarly one defines the rectangle Ql^ the only difference being that the vertical sides 
contain now \{2L + 1)2+^] sites. 

Definition 3.1. A distribution P of b.c. for a rectangle R (which will be Rl, Ql or 
some translation of them) is said to belong to V(R) if its marginal on the union of 
North, East and West borders of R is stochastically dominated by (the marginal of) the 
minus phase tt^ of the infinite system, while the marginal on the South border of R 
dominates the (marginal of the) infinite plus phase tt^. 

The most natural example is to take P concentrated on the boundary condition 
T = — 1 on the North, East and West borders, and r = +1 on the South border. 

Definition 3.2. For any given LgN, 5>0,t>0 consider the Ising model in Rl, with 
boundary condition r chosen from some distribution P. We say that A{L, t, 6) holds if 

E||/if -vr^l <5 (3.1) 

for every P G V(Rl). The statement B(L^t^S) is defined similarly, the only difference 
being that the rectangle Rl is replaced by Ql (and P is required to belong to V{Ql) )• 

With these definitions the iterative scheme developed in [34] can be summarized as 
follows. 

Proposition 3.3 (The starting point). For every /3 (thus not necessarily large) there 
exists c — c{/3) such that for every L G N the statements A{L,t,e ) and 

BiL.t.e-^^'""^'^^^') hold. 

Remark. Notice that the factor g-^^^^^^^ in front of the time t is nothing but the 
negative exponential of the shortest side of the rectangle. 

Theorem 3.4 (The inductive step). For every /3 large enough there exist constants 
ci, C2, C3 such that: 

A{L, t, 6) ^ B(L, ti, 61) ^ A{2L + 1, ts, ^2) (3.2) 
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Q, 



Figure 3. The box Ql and its covering with the rectangles A, B. 

where 

6i = ci (5 + + L^e-^^^^s^) ; h = 2t (3.3) 

^2 = ci(5i + = C3(5 + + L^e-^^^^s^) ; ^2 = e^^^'^h = 2e^3^''t (3.4) 



Remark. In the original statement in [34| the obvious requirement of /3 large was 
missing due to a typo. 



Corollary 3.5 (Solving the recursion). In the same setting of Theorem 3.4, for every 
L G {2^ — IjnGN there exists 



A(L) < exp [-cL'^^j (3.5) 

such that A{L,t,A{L)) holds for every t > T{L) := e^^"'. 

In turn, at the basis of the proof of Theorem |3.4| besides the so called Peres- Winkler 
censoring inequality (see [38j and [34j Section 2.4]), there were two key equilibrium 
estimates on the behavior of (very) low temperature Ising interfaces which we now 
recall and which were the responsible for both the various e~^^ error terms in 61,62 
and the constraint /3 » 1 on the inverse-temperature. The latter was necessary since 
the techniques of ^34| were based on several results of |11| on the Wulff construction 
which in turn use in an essential way low temperature cluster expansion. 

3.1. Equilibrium bounds on low temperature Ising interfaces used in f34|. 

The first estimate is the key for the proof of the first part of the inductive statement 
namely A{L,t,6) ^ B{L,ti,6i). Given the rectangle Ql write it as the union of two 
overlapping rectangles, each of which is a suitable vertical translate of the rectangle Rl 
(see Figure [3]). Call B the lowest rectangle and A the highest one. Then 

Lemma 3.6 (see Claim 3.6 in |34 ). There exists c = c(/3,£:) > such that 

J2 E = +) - 7r-'-(a, = +)) < e"^^'^ (3.6) 

where ti^~ denotes the Gibbs measure in Ql with minus boundary conditions on its 
lowest side and r on the other three sides. 



QUASI-POLYNOMIAL MIXING OF 2D ISING WITH PLUS BOUNDARY 



15 



92 



^ D D 

Figure 4. Two possible arrangements for the open Peierls contours in 
the Gibbs measure under a ( — , +, A) boundary condition. 

In turn, by suitably playing with monotonicity properties of the measure tt as a 
function of the boundary conditions (see the short discussion in the proof of Claim 3.6 
in [34J), the proof of the Lemma can be reduced to establishing the following bound. 

Consider the enlarged rectangle El with sides 3L and 2[(2L + 1)2+^] respectively, 
which can be viewed as consisting of six rectangles Ql stacked together. Let 7r(~'~'+'~) 
be the associated Gibbs measure with (— ) boundary conditions on the North, East and 
West sides and (+) on the South side. For any spin configuration a G { — 1, +1}^^ let 
7 = 7(a) denote the unique open contour corresponding to these boundary conditions. 
Then 

Lemma 3.7. For any /3 large enough there exists c = c{/3^ e) such that for any L 

7r(" (^7 reaches height L^+^) < e""^^'' (3.7) 

Notice that the height L^^^ is well beyond the typical 0{y/L) fluctuations of the 
interface. 

The second equilibrium bound is required for the proof of the statement B{L^ ti. Si) ^ 
^(2L + 1, ^2, S2) (see Section 3.2 and in particular Claim 3.10 in [34j). Here the bottom 
line is the following bound. 

Let Rl consists of two copies of Rl stacked one on top of the other. Let A C ORl 
consists of those boundary sites x = (z, j) in the South border such that \i — L/2\ < ^L^^ 

and j — 0. Consider the Gibbs measure ^^^'^'^^ on Rl with (— ) boundary conditions 
on the union of the North boundary and A and (+) on the rest of ORl- Let Fi be the 
event that the open contour 71 starting on the upper left corner of Rl ends at the left 
end of the interval A without ever crossing the vertical line at z = L/2. Then 

Lemma 3.8. For any /3 large enough there exists c = c(/3, s) such that for any L 

4'^'^^(n)<e-^" (3.8) 

In the scheme envisaged in [34] the role played by the tiny extra piece of (+) boundary 
conditions at the vertices of A, being the main source of the e^^^^^ factor relating the 
time scales ^2,^1 in ( |3.3| ), is quite crucial and therefore it needs a bit of explanation. 

Let us first explain why the length of A was chosen equal to L^^ . Under the boundary 
conditions (— , + , A), for any configuration there exist exactly two open Peierls contours 
71,72 with two possible scenarios for their endpoints (illustrated in Figure [i]): 
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Figure 5. The box i?2L+i and its covering with Qf""^', Q^f\ Q'l^^\ 
bold the exceptional set A. 



In 



(a) 71 joins the two upper corners of Rl and 72 the two ends of the interval A; 

(b) 71 joins the left upper corner of Rl with the left boundary of A whereas 72 joins 
the right upper corner of Rl with the right boundary of A. 

In [34] it was shown, using a significant part of the main machinery of |11 , that the 
ratio between the probabilities of the two cases is roughly of the form e"^^^^!^!"^^)'^/^^^) 
where D is the Euclidean distance between the left upper corner of Rl and the left 
boundary of A. Clearly D ^ L/2 — \A\ ^ 0{Lp'^) and therefore case (b) is much more 
likely than case (a) iff |A| L?^ . The choice L^^ was clearly not optimal and just a 
very safe one. Once the first scenario can be neglected then the fact that 71 does not 
intersect the vertical line at z = L/2 is quite natural (but painful to prove). 

Next we sketchily explain why the need of attracting the contour 71 deep down inside 
the rectangle R^. 

When proving the implication B{L^ti^5i) ^ A{2L + 1,^2,52) we can imagine that 
the rectangle R2L+1 is written as the union of three copies of the rectangle Ql denoted 
by Qf'^^Q^f^QL^^^ (see Figure [sj. 

For simplicity suppose that the boundary conditions around i?2L+i are the "extreme 
ones" namely ( — ,—,+,—) ordered clockwise starting from the North one and imagine 
starting the dynamics from all pluses. 

The Peres- Winkler results allow us to e.g. first run the dynamics in the central 



rectangle for a time ti and then in the left and right ones for some other time 

lag. Thus the dynamics in Q^f^*^ runs with b.c. ( — ,+,+,+) and after a time lag ti 



it will be close to the Gibbs measure vr^ceir^'^^ by less than 61 because of B{L,ti,6i). 

The trouble is that the marginal of this measure on e.g. the East boundary of Q^^* is 
not dominated by tt^ because the unique open contour joining the left upper corner of 



Qcentr ^^le right one will stay close to the upper side of . Therefore we cannot 

use statement B{L^ ti, 61) for the dynamics in Q^^^ to force equilibrium there in another 
time lag ti. 
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An appealing and very intuitive possible way out of this serious problem would be to 
run many times the dynamics in Q^f^*^ until a large deviation forces the open contour to 
go below and to the left of the East side of Q^^*- Since the probability of this fluctuation 
is 0(exp(— cL^^)) it would be enough to wait 0(exp(cL^^)ti) runs. However a rigorous 
implementation of this idea is far from trivial and in ^34j the solution was another one, 
less natural but much easier to carry out. 

If one, by brute force, flips the boundary conditions inside the interval A on the 
South side of Q^f^^^ (— ) the mixing time of the dynamics cannot change by more that 
exp(c(/3)|A|) (see [34| Section 2.5] for more details). Once the boundary conditions have 
been flipped then, thanks to ( |3.8 ), the contours in Q^f^*^ will follow scenario (b) above 
and the resulting distribution over the East boundary of Q^^* will now be dominated by 
the minus phase tt^ allowing another application of the inductive statement B(L^ ti. Si) 
to Q'f' and Q2^^\ 

4. A NEW RECURSIVE SCHEME 

In this section we modify the recursion scheme of |34 and, modulo two equilibrium 
estimates very similar to Lemma 3.7 and 3.8, we prove Theorem [T] We begin by fixing 



some notation. 

Let G N be a large integer, let L = L^y = 2^ — 1 and choose A'o to be the 
smallest integer such that L^^ := 2^° — 1 > [log(L)^J. In our recursion A'o and A^ will 
represent the initial and final scales respectively. To any intermediate scale n G [Nq^ N] 
we associate a length scale = 2^ — 1. We also define the rectangles RmQn to 
have sides (parallel to the coordinate axes) of length (L^, knVL^) and (L^, KNy/Ln+i) 
respectively where = V i^N = 0(log(L)-'^/^) and is a positive constant that later 
will be chosen large enough depending on /3. Thus the very definition of the rectangles 
depends on the final scale. It is worth noticing that ^ /^7v\/i^n+i foi" any n G [Nq, N] . 
Finally, for any n E [A'o,A^], we define the statements A{Ln,t,6) and B{Ln,t,6) as in 
Definition lO 

Having fixed the basic notation our inductive scheme can be formulated as follows. 
We repeat the result on the starting point for completeness, despite it being completely 



obvious after Proposition |3.3| and the remark after it. 

Proposition 4.1 (The starting point). For every /3 there exists c = c(/3) such that for 
every n ^ [A'o,A^] the statements A{Ln^t^ e~^^ ""^^ ^'^) and B{Ln^t^e~^^ ""^^ ^'^) hold. 

Theorem 4.2 (The inductive step). There exist constants ci,C2,C3 and for every /3 > 
/3c there exists such that for any > Hio, for any N large enough and for any 

A{Ln, tn, Sn) B{Ln, 4, 5^) ^ ^(L^+i , t^^+l , 5n+l) (4.1) 

where 

5; = ci(5„ + e-^2'^^+L2e-^^l°g*") ; 4 = 2t„ (4.2) 

^n+i = C3(<5„ + e-^2'^^+L2e-^2i°g*") ; = e^^'.^^tn (4.3) 

Corollary 4.3 (Solving for the final scale). In the same setting of Theorem \4-^ there 
exists c > such that, if t^ '-— e^'^^ and 6]\i := ce~^ '^^ , then for any N G N large 
enough statement A{L]s[^t]sf^ 5 n) holds. 
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Proof of the Corollary. Choose t^^ = ^c' k,n'^ j^^. gQ^^g > 0. Thanks to Propo- 
sition |4.1[ for any /3 > it is possible to choose c' — c '{/3) in such a way that 



^(Ltvq, ^atq, ^A^o) holds with — e ^^^'^l^. Theorem 4.2 immediately implies (use 
> ^ATo) that tfq < e^'^^^ and 6]\i < ce~^ for some other constant c. □ 



Once Corollary 4.3 is proved, Theorem [T] and its corollaries (Corollaries [2] and [4]) 



follow by exactly the same arguments envisaged in [34] for the analogous results. 



In turn the proof of Theorem 4.2 follows step by step the proof of Theorem |3.4| in 



34 once we assume two key bounds on Ising interfaces that we state below. 



4.1. Two key equilibrium estimates for the new recursion. Consider a rectangle 
with boundary conditions that are identically (+) on the South boundary and (— ) else- 
where. The following proposition addresses a large deviation estimate for the vertical 
fluctuations of the unique open contour in this setting, as illustrated in Figure [6) 

Proposition 4.4. Let R be a rectangle of dimensions £ x a^/I with 1 < a < y/I and let 

TT^ ' ^ be the corresponding Ising Gibbs measure with ( — ,—,+,—) ordered clockwise 
starting from the North side. Let A = A(cr) denote the unique open Peierls contour of 
the spin configuration a G ^r. Then for any /3 > /3c there exist constants ci,C2 > 
depending only on /3 such that for any < S < 1 and a, i as above 

^5^"'"'+'"^ (^a : A(a) reaches height 6aVI^ < ^^ie-^2(M' . (4.4) 

Remark. In the proof of the statement ^(L^,t^,5^) ^ B{Ln^t'^^5'^) for n G [A^o,^], 
the above proposition is used with I — 3L^, a — 3^/2kn and 6 = {3^/2)~^. Thus for k 
large enough depending on /3 and for every n G [A^o, N] the r.h.s. of (|4.4|) is quite small. 
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Figure 7. Open Peierls contours confined to the left and right halves 
of the rectangle R under ( — ,+, A) b.c, addressed by the equilibrium 
estimate of Proposition |4.5[ 



The second equilibrium bound that is needed can be formulated as follows. Mark 
the rectangle R as given above by the corners [x^y^y'^x') clockwise starting from the 
Northwest corner. Consider the Ising Gibbs measure on R with the following b.c: 

(i) —1 on the North boundary and on an interval A of length sa^ belonging to the 
South boundary and centered around its midpoint; 

(ii) +1 elsewhere. 

We refer to these boundary conditions as the b.c. (— , +, A) and let u and v denote the 
West and East endpoints of the interval A centered on the South border. 

In this new setting we aim to show that w.h.p. in the random-line representation 
there are two open contours Ai and A2 with 5\i — {x, u} and 8X2 — v} and such 
that Ai (resp. A2) lies entirely in the left (resp. right) half of i?, as shown in Figure [t] 
(recall the discussion following Lemma |3.8| for the role of this event in the inductive 
scheme). This is established by the next proposition. 

Proposition 4.5. For any /3 > /3c there exist ci, C2, > depending only on (3 so that 
the following holds. Let R he a rectangle of size i x ay/i with 1 < a < (l/so)\/I and let 
A = [u,v] be an interval of length sa^ centered on the South border for some s > sq. 
Let V denote the event that there are two open Peierls contours confined to the left and 
right halves of R and connecting the top corners with u^v. Then 
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Figure 8. Infinite strip with mixed b.c. and its unique open contour. 



Remark. In the proof of the statement >B(L^, t^, 5^) ^ ^(L^+i, t^+i, the above 

proposition is invoked with a choice of ^ = L^, n G [A^o,^], and a = 2\/2/^7v, so that 
a = o(\/Z); also, the r.h.s. of (4.5) is always very small provided that the constant Hi is 
chosen to be large enough. 



5. Equilibrium crossing probabilities for the infinite strip 

In this section we study the behavior of the unique open contour in the infinite strip 
with boundary conditions (+) in the upper half-plane and (— ) in the lower half-plane. 
Deriving sharp estimates for the probability that this contour is confined to the upper 
half-plane, as well as a large deviation estimate for the its vertical fluctuations, will 
later serve as a key element in the proofs of Propositions 4.4 and 4^ The analysis 
in this section hinges on the duality tools developed in |40[|41| , which enable us to 
characterize the Ising interfaces for any /3 > /3c. By using this machinery together with 
some additional ideas we establish various properties of the contours, roughly analogous 
to Brownian bridges with logarithmic "decorations" . 

For S dl? define the boundary condition 77 G {±1}^^^*^ to be 



-1 y>0, 
+1 y<0. 



We focus on the case where S is the infinite-strip of width ^, 

5= {1,...,^} X Z, 



(5.1) 



(5.2) 
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whereby the above b.c. 77 gives rise to a unique open contour A connecting the dual 
vertices + h^} (^^^ Figure Such contours have been intensively 

studied and the scaling limit of A is known to be the Id Brownian bridge between 
these two points f20] while our proof requires more quantitative estimates. Tight large 
deviation estimates for vertical fluctuations of A are necessary in several places in our 



proof. This will be established by Theorem 5.3 below (in a slightly more general setting 



via an argument akin to those used for controlling the deviations of the Brownian bridge 
yet carried out within the duality framework of |40|[4T]. 

Significantly more delicate is the crucial estimate of obtaining a lower bound on the 
probability that A is contained in the upper half-plane. The Brownian bridge heuristic 
suggests that this event holds with probability proportional to l/£, and as the following 
theorem confirms this is indeed the case. 



Theorem 5.1. Let S be the infinite strip of width £ with b.c. rj as given in ( |5.1[ );( |5.2[ ). 

For an Ising configuration a on S let \ — A(a) be its unique open contour in the dual 
lattice 5* (i.e. SX = {(|, ^),(^+ |, I)})- For i ^ Z let ^ {\, \} x {i ^ \} 

be the vertices comprising the i-th horizontal level of S"" . Then for every i, 

c C 
- < TT^ (a : A(cr) stays above < — , 

where c, C > are constants that depend only on /3. 

As an immediate consequence we obtain the following lower bound on the spin-spin 
correlation at high temperature for two points on the horizontal boundary of the half- 
strip 5+ = {^,. . . ^} X Z+*, which to our knowledge was previously unknown. 

Corollary 5.2. Let u — {\^\) ciTid v = {£ + ^'^)- every (3 > j3c there exist 

constants c, d such that 

Note that the above corollary also extends to other geometries, for instance rectangles 
with a wide range of aspect ratios where v correspond to the upper corners. We 



postpone the proof of Theorem 5.1 and Corollary |5.2| in order to first obtain several of 
the ingredients that it would require, the first of which being the aforementioned large 
deviation inequality for the open contour in the infinite strip S. 

Theorem 5.3. Let S = S{a,b) be the infinite strip {1, ...,£} x Z excluding the hori- 
zontal slits {1, . . . , a} X {0, 1} and {b, ...,£} x {0, 1} for 0<a<b<£^l with b.c. 77 
as defined in ( |5.1| ). For an Ising configuration a on S let \ — A(a) be its unique open 
contour in the dual 5*; and for i ^ X let = {^,...,^ + ^} x {z+ Then there 
exist some constant C(/3) > such that for any £ the following holds: 

7r| ^cr : A(cr) reaches H^^^ < C exp (— /^^ x^) for all x < y/I ^ 

TT^ ( cr : A(cr) reaches J < C exp {—1^(3 h) for all h > £ ^ 



where > is the constant in the sharp triangle inequality of the surface tension Tf^. 
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Remark. It is fairly straightforward to establish upper bounds as above with an extra 



prefactor of order \b—a\ (see e.g. the first inequality in (5.12)). Eliminating this spurious 
prefactor requires a delicate multi-scale analysis. 

Proof of theorem. In what follows we will prove the following inequality, which is a 
stronger form of the required large deviation estimates: For some C = C{/3) > 0, 



4 



A((t) reaches H^) < Ce ''^yb-a-i 



, h) 



for any h > 



(5.3) 



(we may clearly assume that b > a -\- 1 otherwise the unique open contour is trivial). 
Indeed, the above probability estimate is clearly increasing in the value of 6 — a, which 
in turn is guaranteed to be at most ^ + 1 (reflecting the bounds in the proposition). 
Notice that by choosing C to be appropriately large we need only consider h > \/h — a. 
Fix some large cutoff height n > {h V i)'^ and let 

Sn = Sn{Zx {-n,...,n}) 

be the strip S truncated at ±n with boundary conditions analogous to 77, i.e. negative 
on the upper half-plane and positive elsewhere. Due to the uniqueness of the Gibbs 
measure on 5, the probabilities we seek to bound are obtained as a limit of the corre- 
sponding ones for Sn diS n ^ 00. Further let = (a + ^, ^) and = (6 — ^, |) denote 
the endpoints of the unique open contour in 5*. Define the height of this open contour 
A C 5* at the horizontal coordinate xG{^,...,^+|}tobe 

ht(A, x) = max{^ : (x, y) ^ X} . 

The main effort in the proof will be devoted to the analysis of the vertical fluctuations 
of the contour A within the inner strip with x-coordinates {a + |, . . . , 6 — |}. It is 
the case that large vertical fluctuations in the margins {i.e. large values of ht(A, x) for 
X < a OT X > b) are far more unlikely and can be estimated via standard properties of 
the surface tension. To control the dehcate fluctuations of ht(A,x) for a < x < 6 we 
will apply a multiscale approach, repeatedly bounding the deviations at the horizontal 
midpoints in a nested dyadic partition of the interval between u and v. 

The first step in the proof is to bound the event that the contour includes a given 
point w = (x, h) G 5* in terms of its coordinates h and x. First notice that by (2.8), 



TT^ {a : w e A(cr)) 



'A: SX={u,v} 



I 



J2 is^(^^ 

A: 6X={u,v} 



(5.4) 



Consider the numerator in the last expression: Corollary |2.5| implies that 



A: SX={u,v} 

wex 



and together with Lemma 2.1 we deduce that for some cq = co(/3) > 



Co 



A: SX={u,v} 
weX 



u — w\\v — w\ 



: exp { — ri3{u — w) — rp{v — w)) 



(5.5) 
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To estimate the denominator in ( |5.4[ ) recall Eq. ( |2.7[ ) according to which 



E 

A: (^A^jw,!;} 



As it follows from GKS that decreasing our domain can only decrease the spin-spin 
correlations, letting = {a, . . . , 6} x {— n, . . . , n} (i.e. Sn is the result of "pushing" 
the West and East boundaries of Sn to a and b resp.) we have 

(cTu CTv) > 7T*g* (cFu CFy) . 

By Eq. (2.3) there exists some ci = ci(/3) > such that the spin-spin correlation 
between v in the dual to the infinite strip 5 = {a, . . . , 6} x Z is 

7Ts-{(Ju(Jv) = ^^^^=exp{-Tf3{\u - V\)) , 

- ^1 

where the o(l)-term tends to diS \u — v\ ^ oo. Due to the strong spatial mixing 
properties of the high temperature region /3* < /3c, the value of ttL {au cry) converges to 
the above r.h.s. exponentially fast in n. Already for n > £^ we could absorb the error 
in the constant ci and obtain that for some c[{^) > 0, 



A: 5X={u,v} 



u 



■.exp{-rp{\u - v\)) . 



(5.6) 



By combining ( |5.4| ) with ( |5.5| ) and ( |5.6| ) we conclude that for some C2 = C2(/3) > 0, 

: exp {—T^{u — w) — rp{v — w) -\- t^{u — v)) . 



{a: A(a)) < 



C2\/\U-' 



u 



w\ \v 



w\ 



(5.7) 

At the same time, by the sharp triangle inequality property (2.2) of the surface tension, 
Ti^{u — w) + T^{v — w) > rp{u — v) -\- Hip {\w — v\ + \u — w\ — \u — v\) . (5.8) 
Recalling that w is at height h it is easy to verify that 

4/i2 



\w ■ 



v\ + \u ■ 



w\ 



V > 



u 



Set C 



1+^5 

h2 



> I and now observe that whenever < \u- 



-v\'^ the last expression is at 



least otherwise it is at least ^h. Using this bound for the r.h.s. of Eq. (|5.8|) 

now allows us to produce the following bound out of Eq. (|5.7|): 

TT^ {a : w ^ X(cr)) < — " ' exp 



w\ \v — w\ 



\u—v\ 



A h 



(5.9) 



Straightforward applications of the above bounds will now yield the required bounds 
on the height of A along the margins x < a + ^ and x > 6 — ^ as well as whenever 
6 — a is a uniformly bounded. Indeed, by symmetry we may assume without loss of 
generality that x <a^ \ and note that in this case w = (x^h) satisfies \w — v\ > \u — v\. 
Applying (5.7) combined with the sharp triangle inequality as in Eq. (5.8) we get 

C2 



TT^ (a\we A(cr)) < , 



w\ 



: exp {—Kp\u — w\ 
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Summing the last expression over all w — [x^y) with x < a + ^ and y > h gives that 
7r|^ (cr : ht(A(cr), x) > h for some x < a + ^) < Ci e'^^^ (5.10) 

for some Ci = > 0, and analogously 

7r|^ {a : ht(A(a), x) > /i for some x > b - I) < Ci e'^^^ . (5.11) 

We now turn our attention to the main task of bounding the vertical fluctuations of 
A along the interval (a + ^,6 — ^). First observe that (5.9) immediately provides the 
bound we seek (Eq. ( |5.3[ )) in the special case where \u — v\ = 0(1) (with an implicit 
constant that may depend on /3): In that case a simple union bound over w = (x, h) 
for X G (a + ^, 6 — ^) yields 



7r| (a : ht(A(cr), x) > h for some a+^<x<6 — |) 



< C2\u - v\ exp (-^Kp A /i)) < C2 exp (-k^j A h 



\u—v\ 



(5.12) 



where C2{(3) > incorporates the uniform bound on \u — v\. Combined with (5.10) 
and ( |5.11 ) this concludes the bound in Eq. (5.3) when \u — v\ = 0{1). 

Let M > 2 be some fixed integer whose value will depend only on f3 and will be 
specified later. Justified by the above argument, assume without loss of generality that 



\u — v\ > and exp — ^'1''^^^^ > 



u — v\ . 



(5.13) 



We claim that this in turn narrows our attention to proving Eq. (5.3) for h satisfying 

(5.14) 



M\u-v\<h^ < ^lu-vf/"^. 



To see this recall first that the lower bound on h is justified by selecting a suitably 
large constant C(/3) in Eq. (5.3). For the upper bound, note that if /i^ > — v]^/"^ 

(in which case > — v\^/^ whereas h > — v\^/^) then (5.13) implies that 

— i;| is at most exp A /i)^ and hence Eq. (5.3) follows from a union bound 

over X G (a + ^, 6 — ^) as in (5.12). 

Consider the event whereby the contour A visits a point w ^ Sn given by 

w — (x, y) for some x G (^a + jj\u — v\^b — jj\u — v\^ and y > h . 
Clearly ^y\u — w\\v — w\ > jj\u — v\ thus we can rewrite (5.9) as 



C3 



TT^ {a : w e A(cr)) < 

'^^ y^\u-v 



(-1^/3 A y 
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where C3 > depends only on ^. Summing over all possible values of y > h we now 
obtain that 

. \u—v\ 



T^§Jcr : ht(A(cr),x) > h) < 



C3 



< C3 Yl + 



— 1; 



C3 



6„ ^ 



y=h y>{h V 'i^D 

y>h 



6 

e"5^/3^ 



6 u 



(5.15) 



<C3exp(-|K^(^ A /.)) 

where the constant C3 > depends only on /3. 

We next wish to extend the above bound on ht(A,x) to hold simultaneously for all 
xG(a+^,6— ^)by means of a dyadic partition of the interval between u and v. Set 



K ^ 



and notice that (5.13) ensures that K > 1. Define the following sequence of refinements 
of the interval between u and indexed by /c = 0, . . . , We begin with the trivial 
partition at level 0, 

Zq ^ — Qj ~\- 2 I ^ — b 2 5 

(k) (k)\ 

and refine level k into level /c + 1 by subdividing each subinterval [Zj_-^,Zj ) into M 
equal parts (up to integer rounding): 



z 



(k+i) _ (k) 



for j 0, 



(k+i) 

M(j-l)+i 



1 + 14 (4^-1 + ^''O + ^ for z = 1, . . . , M - 1 and j = 1, . . . , 



Observe that for all admissible j^k we have 

M-^\u-v\-2 < zf^ -zf\ < M-^\u-v\ + 2, 

where the additive terms account for the rounding corrections along the refinements. 
In particular, the expression in the lower bound on the sub-interval lengths satisfies 

M~^\u -v\> M~^\u -v\> ^J\u - v\ > 10 

(as \u — v\ is large enough). Next, define 

hj, = M-^/^/i for = 0, . . . , K , 

and let T^-^^ be the event that the height of the contour at zj^^ does not exceed J2i<k 
rf^ = J a : ht (A(a), ^) < ^ /i^ Uor A: > 1 and 1 < j < 



Recalling (5.15) and rewriting it in terms of T^-^^ and its complement T^"^^ we have that 



hi 
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Exactly the same argument yields that for general fc, 1 < j < and 1 < z < M — 1, 



hi 



Ak) _ Jk) 



A hk 



(5.16) 



To estimate the last expression, observe that hk/{Zj'''^ — Zj^^i) increases with k roughly 
as M^*^/^. More accurately, 

hk 



< 



< 



Jk) _ (k) - M-k\u-v\-2 ~ \u-v\\ M-k\u-v\-2 

Zj Z-_-^ II II II 



— 1^^ + 



h 3k h 
< -M-- r 

4 \u — v\ 



(5.17) 



(where we used the fact that M — v\ > 10) and similarly 



hk 



> 



4 , ^3fc h 
> -MX 



u — v\ 2 5 \u — v\ 



(5.18) 



Our choice of K and the upper bound (5.14) on enable us to derive from (5.17) that 
for all k< K, 



hk 



Ak) _ (k) - 4 



< 5^3i^/4j^|^_^|-3/8< 5 <1 



V2' 



4v^ 



In particular, this identifies the minimizer of the exponent in the r.h.s. of (5.16) and 
implies that 



TT 



s„ [^M{j-i)+i ' S'-i ^^j ) < exp 



hi 

^ ^ Jk) _ (k) 
^3 



Crucially however, the lower bound (5.18) also gives that 

{k) _ (k) - 5 \u-v\ 5 



u — V 



To simplify the notation put p = 1^/3^;^::^ and recall that p > f/^/^M by (5.14), 
hence we may take M sufficiently large so p would also be large. The combination of 
the above inequalities together with a union bound gives 



K-l 



k=l j 



K-l 



-I U UTf = 4. U U {tS'o-:,.. . . Tf-"} 



k=l ij 



K-l 



k=2 



where we used that p > 2 and \/M > logM for any sufficiently large M and it is 
understood that T^-^^ is the full probability space. 
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We have reached level K at which point we wish to examine the remaining points 
altogether. Fix some x E {zj^^^\ Zj^~^'^) and let 

n^ja: ht(A(a),x) <5^/iA. 

i=i ^ 

As estabhshed before hx/izf^^ — z^jl\) < 1 and so 



-I ' ' ^ ^3 exp ( - h,^^,,^^ ■ (5.19) 

On the other hand, by the definition of K we have that > \u — v\^l^ IM (with the 
factor of M due to possible integer rounding in K) and hence 

W_ W 5 \u-v\ 5 iz-i; 



where the last inequality is due to the lower bound on h in (5.14). It now follows 
from (5.13) that 

exp - ^(K)"^ ^(K)^ ^ (- - ^r^^) < 1^ - ^r^ . 



Together with (5.19) this implies that 



\ J \u — v\ \ ^ \u — v\ 

Summing over at most \u — v\ possible choices for x we may now conclude that 



r| (ht(A, x) > for some x E (a + ^, 6 — ^)) < C4 exp f — j^a^/^ | -A , (5.20) 

where = Ylo M''/^ < (l - M-V4)-\ 

Finally, by choosing M to be sufficiently large we can obtain that aj^ < ^ and 

plugging this in ( 5.20[ ) (while recalling that we are in the regime where < h due 

to Eq. ( 5.14| )) concludes the proof of (5.3), as required. □ 



Remark. The truncation argument that was used in the proof of Theorem |5.3| to 
reduce the problem to a finite domain is applicable in our upcoming arguments as well. 
Henceforth, when needed, we will thus work directly in the infinite volume setting to 
simplify the exposition. 



We now introduce the main conceptual element in the proof of Theorem |5.1[ Recall 
our aim is to show that the open contour in the infinite strip S' = {l,...,^}xZ has a 
reasonable probability — namely of order c/£ — of remaining in the upper half-plane 
{i.e. above the dual line {(x,y) : y = — ^})- 

Our approach, based on the Brownian bridge heuristics, is iterative and very much 
based on the intuitive picture in which the open contour really consists of two sim- 
ple lines 71, 72, traveled at constant speed, one starting from the left boundary and 
moving towards the right boundary and viceversa for the second one, meeting in some 
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Figure 9. The open contour in the infinite strip with mixed b.c, pro- 
gressively exposed as two curves 71,72 originating at its endpoints. 

intermediate point. Such a picture, which can be made more precise by progressively 
revealing the contour from the left to right and from right to left (see Figure [9]), allows 

hitting times kind of arguments that we now explain. Let rj^^ (i = 1, 2) be the hitting 
time of either level or level 2^ for the curve 7^. Then, conditioned to the event that 

(i) 

both curves at their respective times Tj have not yet joined and are both at level 2-^, 
by monotonicity and symmetry, with probability at least 1/4 both curves will either 



hit the next level 2-^+^ or join together before hitting level (see Claim 5.8 below for a 
precise formulation). Thus, with probability at least we can force both curves to 
either hit level 2^ or join together before hitting level 0. However, and that explains 
the heuristic bound 1/^, once the curves are at level 2^ ^ then with probability 
bounded away from they will join together without hitting level 0. In other words it 
is enough to force the curves to climb only n — ^ log2 ^ levels in order not to hit level 0. 

The above sketch, however, suppresses a number of technical difficulties such as 
the boundary conditions and the dependence between the two contours. Moreover, 
and contrary to the behavior of the Brownian bridge, the law of the contour A is in 
fact asymmetric w.r.t. the horizontal axis. This follows from our splitting-rule, which 
introduces a vertical bias for the contour: For instance, as illustrated in Figure [T0| 
applying the SE splitting-rule clearly has the open contour move up with probability 
uniformly bounded away from ^. 

To overcome this difficulty we consider the open contours formed by both the SE 
and the SW splitting-rules, 7se and 7sw resp., and examine their union F = 7se U 7sw 
Most importantly the law of their union is symmetric w.r.t. the horizontal axis. We will 
show that F is essentially a "tube" of logarithmic width surrounding 7se, with added 
"decorations" from 7sw which are components of at most logarithmic diameter (and 
similarly if we reverse the roles of 7se, 7sw)- Up to these logarithmic corrections we may 
implement the heuristics of our above sketch, as stated in the following lemmas. Here 
and in what follows we associate with an open contour 7 going from uio v unit speed 
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Figure 10. Asymmetric behavior of the contour due to the global SE 
splitting-rule. 



parametrization 7(t), justifying hitting-time type of events (e.g. "7 hits the vertex y 
prior to hitting etc.). 



Lemma 5.4. Let S = S{a,b) and iJ* (i G Z) be as in Theorem 5.3. For an Ising 
configuration a on S let'^^^{a) and ^swi^-) be the two unique open contours in formed 



by the SE and SW splitting-rules resp., i.e. going from (a + |, ^) to (6 



\,\). There 

exists some C^{P) > so that for any h > 1 the contour ^se (resp. ^sw) hits H'^^_ 



before hitting either 



-C* lO! 



1 1^ 

2' 2> 



with probability at most ^ + /h. 



-C* \ogh 



Proof. We define the gain of a connected subset of dual edges A in the infinite strip S in 
the interval / C Z over distance m, denoted by gii{A^ /, m), to be the maximal difference 
in ^-coordinates between any two points in A whose x-coordinates are contained in / 
and differ by at most m: 

gii{A^ /, m) := max {\y — y\ : \x — x\ < m, (x, y) G A, {x , y) G A, x, x' G /} . 

(5.21) 

We define the gradient of A as its gain over distance 0. The following claim bounds 
the gain of 7 in a neighborhood of a and b\ 

Claim 5.5. Let 7 be the open contour with either SE or SW splitting-rule in the infinite 
strip S — S(a^b) of side-length i defined in Theorem 5.3. Then for any c > there 



exists a constant — Ci{^, c) > such that for all 1 < m < 

7r| ^ gn(7, [a — m^a + m],c log m) > Ci log < Ci/m . 
The analogous statement holds replacing a with b. 

Proof of Claim [jTjf Define / = [a — m, a ^ m] H Z and in what follows take > c. 
Further let = (a + |, |) and = (6 — ^, |) denote the endpoints of the open contour 
7, and let = C'l be some constant to be determined later. Define the set 



S = {((x,y), {x ^y)) G 5^ : x,x' G /, |x — x \ < clogm, \y — y\ > Cj^logm} . 

If gn(7, /, clogm) > Cj^logm then there exist z = (x^y)^z' = ^ 7 such that 

(z, z^) G S. Taking a union bound over ordered pairs of intermediate points z, z'' G S 
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such that z — (x, y) and z' — [x' ^ y') we get that 



TT^ ( gn(7, /, c log m) > log m ) < 



■{z,z')^'E. A: 5\={u,v} 
z.z'^X 



E 

A: 5\={u^v} 



q-S.{X) 



As we have already seen, Eq. (2.3) provides a sharp estimate for the above denominator 



and it remains to consider the numerator. Recall Corollary |2.5| that treated the measure 
of all open contours A in a domain A* that go between two endpoints v as well as an 
intermediate point z, bounding it from above by the product of the spin- spin correla- 
tions TT^* (au cFz) and TT^* {a^ a^). Following essentially the same proof, (4H Lemma 5.4] 
gives an analogous version of this statement for all such contours A going through two 
ordered intermediate points z, z' (that is, 7 connects u to z, thereafter proceeds to z' 
and ends at v) whereby 



E 

A: 5X={u^v} 
z.z'^X 



(5.22) 



Therefore, Lemma [27T] implies that 

C| exp ( - [Tf3{u -Z)+ T(3{Z - Z') + T^{V - z')\) 



A: 6X={u,v} 
z,z'eX 



u 



z\ • \z 



Since ti^{9) > r^(0) (see e.g. [g]) we can bound the last exponent from above by 

^-rpm\u-z\ + \z-z'\ + \v-z'\] ^ ^-rpm^^-v\+c'(\y\/m+\y-y'\)] ^ 

where > is an absolute constant; indeed, the last inequality is justified by the fact 
that \y — y'\ is at least a constant times \x — (recall that > c) and a similar 
statement holds w.r.t. y/m compared to |x — a|. Since summing over x, amounts to 
a factor of O(m^), absorbing an additional 0{m) term from the sum over y/m while 
recalling that \y — y'\ > C\\ogm now implies that 

qs*{\) < Cm~^\u — v\~^^^ ex.p{—rp{u — v)) , 

A: 6X={u,v} 
z.z'^X 

where p can be made arbitrarily large by taking large enough. In conclusion, 

^|(gn(75^5clogm) > logm^ < Cl/m, 
completing the proof. □ 

Claim 5.6. Let 7se cind 7sw be the open contours with the SE and SW s plitting-rule 
resp. in the infinite strip S — S(a^ b) of side-length i defined in Theorem 5.3. Then there 
exists some C2 = C2(/3) > so that for alll < m < i with probability at least l — C2/m 
every connected component o/7sw\7se with zero distance from ([a — m^a-\-m] xZ)n7sE 
has diameter at most C2 logm. 
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Proof of Claim \5.6[ Let B denote the event that there exists a connected component 
of 7sw \ 7sE with zero distance from ([a — m^a + m] x Z) H 7se and has diameter at least 
logm. We begin by conditioning on 7se. The contour partitions S into two sets 
Stop and S^bot- For a set of dual edges A let V{A) C denote the set of vertices at 
distance ^ from A. The effect of conditioning on 7se is equivalent to conditioning that 
o'u = Vu wh ere U = V{A{js^)) (recall the definition of the edge-boundary A(-) in ^ 



after Lemma 2.2) and rju ^ +1}^ is the configuration given by 



1 U ^ S'top? 

+1 



Conditional on 7se the configuration a on S \ U is given by the Ising model on S^top \ U 
and S^BOT \ U with minus and plus boundary conditions respectively. Let 9 denote the 
ensemble of contours of this configuration given by the SW (not SE!) splitting rule. 
Since the boundary conditions are all minus and all plus there are no open contours. 
Every maximal connected segment of 7sw \ 7se must be a subset of one of the closed 
contours of 9 and must share a common vertex with 7se. 



By Theorem [5^ we see that 7se U 7sw C A* where A = {1, . . . x {-£, ...,£} 
except with probability 0(exp(— c^)) for some c(/3) > 0. By Claim [5^ there are at 
most 2Cim\ogm vertices in X = ([a — m^a + m] x Z) H 7se except with probability 
Ci/m. For z, G A* the probability that both lie in the same closed contour of 9 is at 
most exp{—T^{z — z')) by Lemma [Z6j Since the surface tension achieves its minimum 
on the sphere at t^(0) > combining the above estimates we have that 

TvliB) < Cexp{-c£) + C^/m + J2 J2 exp{-\z - z'\r^{0)) 

zeX z'\\z-z'\>C2\ogm 

< Cexp(-c^) + Cl/m + Crn^-^^^^^^^^ log^ m , 
The desired result follows from a sufficiently large choice of . □ 

Finally we show that the contours 7se and 7sw are unlikely to travel much farther 
than in the x-coordinate before attaining height h or —h. 

Lemma 5.7. Let S — S(a, b). For any w > and < h < £ define the rectangle 

7^={a— | — It;, ...,a— ^+it;}x{^ — /i, 

Let B denote the event that the contour 7se (resp. jsw) beginning at (a — ^, ^) exits TZ 
and the first point it hits inTZ^ is in {a — ^ — w — 1^ a — ^ + w + 1} x — . . . ^ h -\- 
There exists a constant C^{P) > independent of w and h such that 

nl{B)<Clexp{-w/iClh')). 

Proof As before let u = (a + ^, ^) and = (6 — ^, ^) denote the endpoints of the 
open contour 7. Let Br and Bi denote the events that the contour exits to the right 
and left of TZ respectively so that B = Bi U Br. We will examine the case that the 
contour exits to the right and the left case will follow similarly. Fix c(/3) = (SeC/^)^ 



where Cf^ is the constant in Lemma 2.1 For large enough C3 the bound holds trivially 



when w < 2ch^ so assume that w > 2ch^ . We may define M{w,h) G Z such that 



32 



EYAL LUBETZKY, FABIO MARTINELLI, ALLAN SLY, AND FABIO LUCIO TONINELLI 



0<w/{4:ch^) <M{w,h) <w/{ch^) and \cMh^ - {h - a)\ > ^|6-a|. Define the set of 
sequences of points 

S = {(xi, yi), . . . , (xm, Vm) ^S^\\/i, = a + ^ + ich^ , \yi\ < h} 

and note tliat if the contour exits TZ to the right then it must pass from u through a 
sequence of points in S in order and then to b. For a sequence ^ G S we say that A is 
^-admissible if it passes through the points in ^ in order and then returns to b. Note 
also that by the construction of M we have that \xm — {b — ^)\ > ^\b — a\. Taking a 
union bound over sequences in S we get that 



4(e«) < 



E E 

X:6X={u,v} 
^—admissible 



I 



X: 6X={u,v} 



In analogy to Corollary 2.5 and (5.22) one has 



M+l 



gA*(A) < Yl ^A* , 



A: 6X={u,v} 
^—admissible 



1=1 



where we denote u = (xo,yo) and v = {xm^i^Vm^i)- Therefore, Lemma 2.1 implies 
that 



C^+i exp ( - E.=|' Mixi, y,) - 



A: 6X={u,v} 
^—admissible 



'UZl^ li^i^Vi) - {Xi-l,yi-l] 

exp ( — M — r/3(« — v)) 
(3hrJl\b-a\ 



bummmg over the {2h + 1)^ elements of S we have that 



E 

X:SX={u,v} 
^—admissible 



exp ( — M — T/3(i/, — v)) 



and it follows (recall Eq. 2.3) that there exists a constant C{/3) > such that 

ttI {BR)<Ce-^ 



and a similar estimate holds for 7r| {Bl)- Recalling that M > w/{Ach?) and taking a 
suitably large now completes the proof. □ 

We now complete the proof of Lemma |5.4[ Recalling that F = 7se U 7sw we define 



1 1^ 

2' 2> 



to (b- 



1 1^ 

2' 2 



. First observe that 



Ftop to be the highest path in F connecting (a + 
Ftop is indeed well defined. The collection of dual edges F partitions S into two infinite 
components and possibly a number of finite components. To construct Ftop, view the 
upper infinite component as a subset of by drawing a unit square centered at each 
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of its points. Then Ttop is its "horizontal" boundary connecting (a+ ^, ^) to (6 — ^, 
Define Fbot similarly as the lowest path in F. Then 

(Ftop hits before hitting or (6 - ^, ^)) 

<7Ts {^BOT hits before hitting or {b - |, ^)) 

= (Ftop hits i7* before hitting or (6 - ^, ^)) , (5.23) 

where the inequality follows by the fact that Ftof lies above Fbot while the equality is 
by the symmetry of F. It follows that 

{^Top hits before hitting or (6 - ^, ^)) < ^ . 



Lemma 5.7 guarantees that except with probability 0(exp(— /i^)) both contours 7se 
and 7sw hit either H'^^_Q^Yogh -^h-C^iogh before traveling distance order in the 
horizontal direction, and as such we only need to consider the interval [a — /i^, a + /i^]. 
Now set 

C*-8Cr(/3,C2*) + 4C|(/3) 
where Cj^,C2 are the constants from Claim [SlS] and Claim [5^ This guarantees that 

gn (7sE, [a-h\a + h% log h) < 4C^) log h , 

gn (7sw, [a -h\a + h%4C^ log h) < 4C^(/3, 4C^) log h , 
max d ((x, y), 7se) < log h 

x^[a—h^,a-\-h'^] 

and similarly around b with probability at least 1 — C~^/h. In particular, given the 
above event we have that the vertical distance between Ftop and 7se does not exceed 
C~^loghm the intervals [a — /i^, a + /i^] and [b — h^^b+h^] which together with Eq. (5.23) 



completes the proof. □ 

Proof of Theorem lower bound. The proof proceeds by progressively reveal- 
ing the contour 7se. Let wq denote some large constant and let Wi = 2wi-i—2C~^ log Wi-i 



where is the constant from Lemma 5.4, Taking wq sufficiently large it is easily con- 



firmed that Wj > c2^ for some constant c > 0. 

Starting from the left at (^, ^) for j > let Aij be the event that 7se hits H^. before 
hitting or reaching + ^, Similarly starting from the right at + ^, |) let 
Arj be the event that the contour hits H^^ before hitting H^^ or reaching Let 
Aj = Aij nArj and let Bj be the event that the contour 7se hits neither H^^ nor H^.. 

We begin by giving a crude lower bound on the probability of Aq. Let U be the event 
that the spin configuration takes the value (+) for all the vertices (1, j) and (£, j) for 
< j < it;o + 1- This occurs with probability at least 

On the event U the contour 7se(c'") directly passes from (^, |) to (^, it;o + |) and from 
+ ^, ^) to + ^, it;o + It follows that U C Aq and hence, 

7Ts{Ao)>{h-'^f'''^\ (5-24) 



<2 

SO in particular Aq occurs with constant probability 
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We will establish the following claim. 
Claim 5.8. There exists a constant c > such that for all j > we have that 

TTsiAj^i U Bj+i I Aj)>^- c2-K 

If we assume the claim then 

> {\-c2-^) [7Ts{A,) + ns{B,)] 

where the first inequality follows from the fact that Aj and Bj are disjoint. Hence by 
induction and equation (5.24) we have that for any fixed positive integer K 

^(^K+|iog,.) + PiBK^^o,,i) > P{Ao) n (i - ^2-^') > c^4-^^-^ . (5.25) 



We now prove the Claim 5.8, Note that if a' > cr then the curve 7se(c'"0 must lie on or 
above 7se(c'") and hence the event Aj^i is increasing in a and so is Aj^iUBj^i. Through 



a series of monotonicity arguments we will relate this event to that in Lemma 5.4 



Suppose that Aj holds and that the left part of 7se first hits H^. at dual vertex 
{zj" + ^, + ^) and denote this part of the contour by 7^". Similarly denote the right 
part of the contour as 7^^ from (^+ |, ^) to {zj^ + ^i^j + |) with ^ < zj" < Zj^ < £ + ^. 
Finally, let Dj^i denote the event that the contour 7se running between (zj" + 1, + ^) 
and {zj^ + ^, + |) either hits H^.^^ at both ends before hitting H^^ or hits neither 
i7 * . , , nor i7 * 1 . With these definitions we claim that 

TVS {Aj+1 U Bj+i I Aj) > Tv^s {Dj+i I au^ = r]u^) , (5.26) 

where tt^ denote the measure on the strip S with boundary condition given by (+) up 
to Wj and (— ) above wj, 

Uj = {[l,zf + 1] U [z^J]) X {wj,Wj + 1} 

and 

{ + 1 U2 Wj, 
1 I 1 

-1 U2 = Wj + 1 . 

The sequence of monotonicity arguments establishing ( |5.26[ ) is best explained schemat- 
ically, see Fig. jJ_ and its caption. 

Let Dj"^^ denote the event that the contour 7se from (zj" + + ^) (resp. (zj^ + 
^, Wj + ^)) hits H^.^^ or {zf + |, it;^ + ^) (resp. {zf + l,Wj + ^)) before hitting H"^^. 
Clearly = ^j+i ^ ^j+i symmetry we have that 

Since these events are both monotone, by the FKG inequality we have that 
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Figure 1 1 . Monotonicity transformations to reduce the two segments 



of the open contour in the infinite strip to the setting of Lemma |5.4 
First (top-right) impose extra (— ) spins at height wj + 1, making the 
increasing event ^j+i U Sj+i more unhkely. Next (bottom-right) set 
the spins at height wj and above the contours if^lf to (+); this does 
not change the probabihty of ^j+i U Bjj^i due to the Markov property 
of the Gibbs measure. Finahy (bottom-left) remove the constraint that 
the spins just below 7^", 7^^ are (+), again making the event Aj+i USj+i 
more unlikely. The dotted line is at height Wj + 1/2. 

It finally suffices to note that the conditional event D^.^ is exactly the event considered 



in Lemma 



5.4 



once we shift the strip down by Wj and set h — Wj — C^^logWj, a 
and b = — 1. Then we have that 

4 (Dj^+i I au, =rju,)>\- 0{l/w,) > ^ - c2-i , 
which completes the proof of Claim |5.8[ 
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To complete the proof of the lower bound in Theorem 5.1 we will show that 

^(^K+l+|l0g2^ I ^X+|l0g2^) > 2 ' 

for a large enough constant K — K{(3). Applying the same monotonicity transfor- 
mations as in Claim |5.8| this reduces to the probability of a contour reaching height 



Wj-^i—Wj > c2^\/Z, which is less that ^ for large enough K by Theorem 5.3 Combining 
with equation ( |5.25[ ) it follows that 

as required. □ 

Proof of Theorem upper bound. Let M. = {w = {x,y) : y > 0} denote the 
upper half-plane and consider the correlation between the spins at u = (^5 |) 



<2' 2> 

V = (£+1, 1) in H*. It is known (see e.g. (35} pl61 Eq. (5.29)]) that for some Ci(/3) > 0, 

By GKS we can reduce the domain to the dual half-strip of width |ix — z;] 
S+ = M*n{w = {x,y) : ^<x<i+^}, 

and obtain that 

vth* (cTu cTv) > vr^+ (cFu (^v) = - ' 

6j={u,v} 6j={u,v} 

where S = {w = (x^y) G Z2 : ^ < x < ^ + ^} is the dual strip and the last inequality 
is justified by Lemma [Z2} Combining with Eq. (|2.3|) we can conclude that 



ttI. (a : Xia) stays above H* ^ = ^'^^=^"'"^ , , < ^\ , ^ = C/i , 

thus completing the proof. □ 

Proof of Corollary | 5.2[ The upper bound follows directly from the GKS inequalities 
and the exact solution |35| in the infinite half-plane. As for the lower bound, one has 

XeS+:SX={u,v} XeS+:dX={u,v} 

^xeS+:SX=^uM 



where in the first inequality we applied Lemma 2.2 and in the second one we used (2.7) 



and the lower bound of Theorem 5.1 Recalling (2.3), the desired result follows. □ 
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6. Proofs of Propositions \4M and [4^ 



The remaining part of this work is devoted to the proof of the two equihbrium 
estimates needed for the new recursive scheme (detailed in SectionQ using the estimates 
obtained thus far. 



6.1. Proof of Proposition |4.4| . Let u,v be the initial and final point of A and as 
before let denote level h of the rectangle i?*. By Eq. ( |2.8| ) we have the following: 

^cr : A((j) reaches height SaVi^ 

J2 iR^w/ E (6-1) 



A reaches if* ^ 
(5A={^i,^'} 



Observe that by monotonicity that (6.1 ) is increasing in the height of R and so without 
loss of generality we take 

i? = 5+ = {(xi,X2) : 1 < xi < ^ and X2 > 1}. 
We now bound separately the numerator and the denominator in (|6.1|). By the same 



argument used to prove equation ( |5.9| ) (note that here we bound the probability of the 
contour exceeding height 5a\fl < ^ by the assumptions on 5 and a), there exists some 
ci(/3) > such that for any z G H^^^, 

J2 9fi'(A) < yexp(-r^(0)^-K/3(5a)2), 

A: 6X={u,v} 



and therefore 



E E 

zeX 



qn-W < ci exp ( - t/3(0)^ - t^/siSaf) . 



(6.2) 



Next we bound from below the denominator in (6.1). Recall that S is the infinite strip 
S = {{xi,X2) : l<xi<i}. Then by Lemmal2?2p 



E 

xcK 

6X={u,v} 



Yl Aci^* qs-W 

SX={u,v} 

Yl Ac^* qs-W 

6X={u,v} 



E 

ACS'* 
6X={u,v} 



qs- (A) . 



(6.3) 



The last factor is estimated in Eq. |2.3| while the first factor can be interpreted as the 
probability in the canonical ensemble given by the weights qs*{') and conditioned to 
start at u and to end at v that the contour A stays above the line at height —1/2. By 
Theorem 5.1, this is of order c/l. The desired claim then immediately follows. □ 



6.2. Proof of Proposition 



4.5 



Let R be the I x a\/J rectangle given by the endpoints 
(x^y^y^x^) clockwise starting from the Northwest corner, and let u and v denote the 
West and East endpoints of the interval A centered on the South border, as was shown 
in Figure [t] in ^ Recall that in our setting we have a b.c. r] which is (— ) in the North 
and on A and (+) otherwise, and that the event V states that A is connected to the 
North via two contours confined to the left and right halves of R respectively. 
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Figure 12. Domain decrease from the rectangle i? to a disjoint union 
of rectangles separating the two open contours. 

Our first step in establishing that V occurs except with probability ci£^^ exp(— C2«^) 
is eliminating (except with the aforementioned error probability) the scenario where 
open contours connect x to y and u to v. 

Lemma 6.1. Let R = {x,y,y\x^) and A be an interval of length sa'^ centered on the 
South border x'y' . For any /3 > /3c there exist C3, C4 > and > depending only on 
/3 such that if s > sq then for any £ we have 

= {x, y} , 6X2 — {u, v}) < exp(— C4a^) . 

Proof. Our starting point is the equality ( |2.8[ ), which allows us to rewrite the probability 
of certain contours in terms of their weights: 



where 



E 



{x,y} , 6X2 = {u,v}) 
QR^ (A) , 



^^2 := 



A=AiUA2 
6Xi={x,y}, SX2={u,v} 



E 

A=AiUA2 
6X={x,y,u,v} 



qr- (A) . 



(6.4) 



By Lemma we have 



E 

Ai 

6Xi={x,y} 



E 

A2 

6X2={u,v} 



where the last equality is by \2.1\ . Plugging in Lemma |2.l| it now follows that 

*i < exp (- [t^{x -y)^ T^{u - v)]) = exp (- [£ + sa^] t^(0)) . (6.5) 

Next we consider ^^2- As before, ( |2.7[ ), the GKS-inequalities ([l8j[23|) and symmetry 
imply that 

^2 = TT^* {cTxCTyCruCrv) > ^R* {crxCru)7TR* {^y^v) = ^/^* {(^x^^uf • 

Furthermore, by GKS, spin-spin correlations are non-decreasing in the domain so 
we can clearly confine our domain to the disjoint union of the two rectangles Q — 
(x, u ^ X ) as shown in Figure 12 and obtain by symmetry that 

^2 > X'^ where X := 7Vg*{axCru) • 
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To control the value of X, let z be the center of the rectangle G and further define 
Qi = (x,a,6, x') and Q2 = (a^u^u^b) to be the left and right halves of each of 



dimensions \{i — sa'^) x aVI. See Figure 13 for an illustration. The GKS-inequalities 
(together with a reduction of the domain) yield 

Another application of (|2.7[) gives 



XcGi XCG2 ^ XcGi ' 

with the equality due to symmetry. Define S to be the infinite half-strip of width 
\{i — so?^ obtained by extending the South border of Q\ {i.e. the edge bx') to — oc. 



Since Qi is a subgraph of S it follows from Lemma |2.2| that 

X>X^ where X := E W • 



We now claim that 



where 



6X={x,z} 



j>l- W^y^ exp(-C2a') , (6.6) 



ACS' 
6X={x,z} 



Indeed, X /Y is precisely the probability that the contour X C Sn whose endpoints 
are 6X = {x,z} stays above the horizontal line x^b (the South border of Gi)- If is 
the midpoint of x and x^ then by monotonicity if we condition all (+) in the rectangle 
{x,a,z,z^) then this only increases the probability that the contour hits the line x^b. 



Proposition|4.4 then establishes equation (6.6). Taking C3 large enough in the statement 



of Lemma |6.1| we can assume that a is at least a large constant times \/log^, in which 
case (6.6) gives X = (1 — o{l))Y. 

On the other hand, taking S to be the strip obtained by extending the North and 
South boundaries of Gi to ±oc, we can now express Y in terms of 

E ^^*(^)- 

ACS' 
6X={x,z} 

To this end, observe that Y/Y is the probability that the open contour A connecting 
x^z in S stays below the horizontal line xa (the North border of Qi). Therefore, by 
monotonicity and Theorem |5. 1| there exists some c = c{l3) such that 

y>c/i, (6.7) 
and on the other hand ^TZ, Formula (2.22)] gives for some > 

F = ^^±^exp(-r^(x-z)). (6.8) 
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Figure 13. Comparing the weight of contours from x to u to that of 
contours connecting both points to the center of the rectangle z. 



Combining ( |6.6[ ),( |6.7[ ) and ( |6.8[ ) it now follows that 

X > d-^i-"^^ exp(-T^(x - z)) , 

and recalling that ^ ^ we deduce that there exists some c''(/3) > such 
that 

^^2 > d'i-"^' exp(-2T^(x - u)) . (6.9) 



To conclude the proof, we combine ( |6.5| ), ( |6.9[ ) and get that 

1^ < (l/c")^^" exp (- {i + sa'\ r;j(0) + 2Tf,{x - u)) . 

Recall that rp{0) is an analytic and even function of for any /3 > /3c, hence in particular 
there exists some c = c(/3) > such that 

r^{0) - Tf3{0) < c9^ for any 6> E R. 
Since in our case 9 < arctan {^0) ^ 2a/\/i it follows that for some d = c'(/3) > 

r^{x — u) < T^(0)|x — u\ + cc? . 

On the other hand, 

i + sc? — \x — u\ — Oi^sc?^ . 
Combining these inequalities completes the proof. □ 

We next wish to show that whenever 8\\ — {x, u) , 8X2 — {y, v} we also have that 
the corresponding open contours are confined to the Ri and i?^, the left and right 
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halves of R respectively, except with an appropriate exponentially small probability. 
The complement event we wish to analyze is illustrated in Figure [M} 



Lemma 6.2. Let R — {x^y^y'^x') and A he an interval of length sa^ centered on the 
South border x'y' . Denote by Ri and Rr the left and right halves of R respectively. For 
any /3 > /3c there exist c^^cq > and sq > depending only on /3 such that if s > sq 
then for any i we have 

TT^ (Ai c Ri, X2 C Rr \ SXi = {x, u} , 6X2 = {y, v}) >1- exp(-C6Q^^) . 

Proof. Let / denote the central column of i?, i.e. I is the vertical line connecting the 
centers of the North and South boundaries of R. We aim to bound the probability that 
the contour connecting x to u crosses /, and similarly for the contour connecting y to 
V. By equality (2.8) we can write the former probability as $i/$2 where 

$1 := QR^ (A) , $2 := 5^ QR^ (A) , 

A=(Ai,A2) A=(Ai,A2) 
(5Ai={x,w} , Ain//0 SX^={x,u} 
6X2={y,v} Sh2={yM 

and a union bound (together with symmetry) gives 

4 (Ai cRi,X2CRr\ 6X1 = {x,y} , 6X2 = {u,v}) > 1 - 2^^ . 

^2 
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To bound $2 from below we compare it to ^^i, defined in (6.4). Indeed, 
$2 



1 



^2 



1 - 7r^(5Ai = {x, y} , 6X2 = {u, v}) >1 - i""^ exp(-C4a^) , 



where the last inequality is precisely the statement of Lemma [671] Combining this with 
the estimate on ^^2 given in (|6.9|) we conclude that for some absolute c'^ > 0, 



$2 > {l-t^exp{-C4a^))c''l 



-c" ^-Tp{x-u)-Tp{y-v) 



(6.10) 



(we can assume that a is at least a large constant times y^log^, otherwise the statement 
of Lemma 6.2 trivially holds). 



Remark. The estimate given in (6.9) for ^^2 in the proof of Lemma |6.1| in fact had an 



exponent of —2ti^{x — u). The above bound featuring —{Tj^{x — u) ^ Tj^{y — v)) in the 
exponent readily follows from (6.9) by the symmetry between {x^u} and {y^v} in the 
definition of ^^2- 

It remains to bound $1. To this end, for a given contour 7 define to be the 
graph with the edge-set i?* \ A(7), where A(7) is the edge-boundary of the contour 7. 
Crucially, in our case the edge-boundary of 7 is disjoint to the edges of Ai since these 
contours are compatible. Therefore, by Lemma [2. 3| 



and plugging this into the definition of $1 we deduce that 



$1 



J2 

7 : 6j={y,v} A : SX={: 



An//0 



The sum over the weights qg^ can be estimated via Corollary 2.5 using a simple union 
bound (see e.g. |8, Eq.(3.3)] for a similar argument). Indeed, 

qg^ W <Y1 - Y ^ii • 



A : 5\={x,u} 



zel A: SX={x,u} 

zex 



At the same time, we can use (2.7) and then increase the appropriate domains to the 



infinite-volume lattice to obtain that 

j:S-f={y,v} 
A : 5X={x^u} 



zel 



An//0 



where here tt^ denotes the (unique) Gibbs measure on at Altogether, 

^1 < ^7r^(cr^cr^)7r^(cr^cr^)7r;^(cr^a^) 



zei 



< Y ^^P (-^/^(^ -x) - T^{z -u) - T^{y - v)) 



zei 
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At the same time, due to the sharp triangle inequahty property of the surface tension 
(Eq. ( |2.2[ )), for any /3 > /3c there exists some a^^^ > such that for any x, z 

Tj^^Z — x) + Ti^{u — z) > Tj^{u — x) + /^^ {\z — x\ \u — z\ — \u — x\) . 

Combining the last two displays implies that 

$1 < e-^M^-^)-^M2/-^) ^ exp [-n^ {\z -x\ + \u-z\-\u- x\)] . (6.11) 



zei 



The first term above cancels when combining (6.11 ), (6. 10) and we obtain that for some 
> 0, 

|i < cV^exp[-/^;3(|z-x| + \u- z\-\u-x\)] . (6.12) 

A straightforward manipulation of the above exponent will now complete the proof: 
One can easily infer from the triangle-inequality and symmetry that 

min (Iz — x\ + \u — z\ — \u — x\) — \v — x\ — \u — x\ , 
zei 

and recalling that the dimensions of R are a VI x i and \u — v\ = sa^^ we have 
\v - = |(^ + sa'^f ^a^l, \u-x\^ sa^f + a^l . 



Hence, as long as a < {1/ s)\/J (guaranteed by the assumptions of Proposition 4.5) we 
have \u — x\^\v — x\ < 2£ and 

sa'^i 1 

n — ^ r > 

x\ + \v — x\ 2 



■ x\ — \u — x\ = -— r > -sa^ . 



Plugging this in (6.12) while summing over the |/| values that z can assume gives that 

<cV+iexp(-(/^^/2)5a2) 



$2 

(with room to spare), as required. □ 

Acknowledgments 

We are very grateful to Yvan Velenik for valuable suggestions on the random-line 
representation and related equilibrium estimates, as well as for pointing out the novelty 
of Corollary |5.2[ EL and AS are grateful to the Mathematics Department of the 
University of Roma Tre for hospitality and support. FM and FLT are similarly grateful 
to the Theory Group of Microsoft Research. FLT acknowledges partial support by ANR 
grant LHMSHE. 

References 

[1] D. B. Abraham and P. Reed, Phase Separation in the Two- Dimensional Ising Ferromagnet, Phys. 
Rev. Lett. 33 (1974), no. 6, 377-379. 

[2] K. S. Alexander, The spectral gap of the 2-D stochastic Ising model with nearly single-spin boundary 
conditions, J. Statist. Phys. 104 (2001), no. 1-2, 59-87. 

[3] K. S. Alexander and N. Yoshida, The spectral gap of the 2-D stochastic Ising model with mixed 
boundary conditions, J. Statist. Phys. 104 (2001), no. 1-2, 89-109. 

[4] A. Bianchi, Glauber dynamics on nonamenable graphs: boundary conditions and mixing time, Elec- 
tronic Journal of Probability 13 (2008), 1980-2012. 



44 



EYAL LUBETZKY, FABIO MARTINELLI, ALLAN SLY, AND FABIO LUCIO TONINELLI 



[5] T. Bodineau and F. Martinelli, Some new results on the kinetic Ising model in a pure phase, J. 

Statist. Phys. 109 (2002), no. 1-2, 207-235. 
[6] J. Bricmont, A. El Mellouki, and J. Frohlich, Random surfaces in statistical mechanics: roughening, 

rounding, wetting..., J. Statist. Phys. 29 (1982), 193-203. 
[7] P. Caputo, F. Martinelli, F. Simenhaus, and F. L. Toninelli, ''Zero" temperature stochastic 3D 

Ising model and dimer covering fluctuations: a first step towards interface mean curvature motion, 

Comm. Pure Appl. Math., to appear. Available at arXiv : 1007 . 3599 (2010). 
[8] F. Cesi, G. Guadagni, F. Martinelli, and R. H. Schonmann, On the two-dimensional stochastic 

Ising model in the phase coexistence region near the critical point, J. Statist. Phys. 85 (1996), 

no. 1-2, 55-102. 

[9] L. Chayes, R. H. Schonmann, and G. Swindle, Lifshitz' law for the volume of a two-dimensional 
droplet at zero temperature, J. Statist. Phys. 79 (1995), no. 5-6, 821-831. 

[10] R. Dobrushin and O. Hryniv, Fluctuations of the phase boundary in the 2D Ising ferromagnet. 
Comm. Math. Phys. 189 (1997), no. 2, 395-445. 

[11] R. Dobrushin, R. Kotecky, and S. Shlosman, Wulff construction. A global shape from local in- 
teraction. Translations of Mathematical Monographs, vol. 104, American Mathematical Society, 
Providence, RI, 1992. Translated from Russian by the authors. 

[12] D. S. Fisher and D. A. Huse, Dynamics of droplet fluctuations in pure and random Ising systems, 
Phys. Rev. B 35 (1987), no. 13, 6841-6846. 

[13] L. R. Pontes, R. H. Schonmann, and V. Sidoravicius, Stretched exponential fixation in stochastic 
Ising models at zero temperature. Comm. Math. Phys. 228 (2002), no. 3, 495-518. 

[14] C. M. Fortuin, P. W. Kasteleyn, and J. Ginibre, Correlation inequalities on some partially ordered 
sets, Comm. Math. Phys. 22 (1971), 89-103. 

[15] G. Gallavotti, The phase separation line in the two-dimensional Ising model, Comm. Math. Phys. 
27 (1972), 103-136. 

[16] R. J. Glauber, Time- dependent statistics of the Ising model, J. Mathematical Phys. 4 (1963), 
294-307. 

[17] L. Greenberg and D. loffe. On an invariance principle for phase separation lines, Ann. Inst. H. 

Poincare Probab. Statist. 41 (2005), no. 5, 871-885. 
[18] R. B. Griffiths, Correlations in Ising ferromagnets. III. A mean-field bound for binary correlations. 

Comm. Math. Phys. 6 (1967), no. 2, 121-127. 
[19] Y. Higuchi, On some limit theorems related to the phase separation line in the two-dimensional 

Ising model, Z. Wahrsch. Verw. Gebiete 50 (1979), no. 3, 287-315. 
[20] O. Hryniv, On local behaviour of the phase separation line in the 2D Ising model, Probability 

Theory and Related Fields 110 (1998), 91-107. 
[21] D. loffe, Exact large deviation bounds up to Tc for the Ising model in two dimensions. Probability 

Theory Related Fields 102 (1995), no. 3, 313-330. 
[22] D. loffe, Large deviations for the 2D Ising model: a lower bound without cluster expansions, J. 

Statist. Phys. 74 (1994), no. 1-2, 411-432. 
[23] D. G. Kelly and S. Sherman, General Griffiths' Inequalities on Correlations in Ising Ferromagnets, 

J. Math. Phys. 9 (1968), no. 3, 466-484. 
[24] D. A. Levin, Y. Peres, and E. Wilmer, Markov Chains and Mixing Times, American Mathematical 

Society, 2008. 

[25] I. M. Lifshitz, Kinetics of ordering during second-order phase transitions, Sov. Phys. JETP 15 
(1962), 939. 

[26] T. M. Liggett, Interacting particle systems. Classics in Mathematics, Springer- Verlag, Berlin, 2005. 

Reprint of the 1985 original. 
[27] E. Lubetzky and A. Sly, Critical Ising on the square lattice mixes in polynomial time, preprint. 

Available at arXiv: 1001 . 1613 (2010). 
[28] F. Martinelli, Lectures on Glauber dynamics for discrete spin models. Lectures on probability 

theory and statistics (Saint-Flour, 1997), Lecture Notes in Math., vol. 1717, Springer, Berlin, 

1999, pp. 93-191. 

[29] F. Martinelli, On the two-dimensional dynamical Ising model in the phase coexistence region, J. 
Statist. Phys. 76 (1994), no. 5-6, 1179-1246. 



QUASI-POLYNOMIAL MIXING OF 2D ISING WITH PLUS BOUNDARY 



45 



[30] F. Martinelli and E. Olivieri, Approach to equilibrium of Glauber dynamics in the one phase region. 

I. The attractive case, Comm. Math. Phys. 161 (1994), no. 3, 447-486. 

[31] F. Martinelli and E. Olivieri, Approach to equilibrium of Glauber dynamics in the one phase region. 

II. The general case, Comm. Math. Phys. 161 (1994), no. 3, 487-514. 

[32] F. Martinelli and A. Sinclair, Mixing time for the solid- on- solid model, Proceedings of the 41st 

annual ACM symposium on Theory of computing (STOC 2009), pp. 571-580. 
[33] F. Martinelli, A. Sinclair, and D. Weitz, Glauber dynamics on trees: boundary conditions and 

mixing time, Comm. Math. Phys. 250 (2004), no. 2, 301-334. 
[34] F. Martinelli and F. L. Toninelli, On the mixing time of the 2D stochastic Ising model with "plus'' 

boundary conditions at low temperature. Comm. Math. Phys. 296 (2010), no. 1, 175-213. 
[35] B. M. McCoy and T. T. Wu, The two-dimensional Ising model. Harvard University Press, Cam- 
bridge, Massachusetts, 1973. 
[36] A. T. Ogielski, Dynamics of fluctuations in the ordered phase of kinetic Ising models, Phys. Rev. 

B 36 (1987), no. 13, 7315-7318. 
[37] L. Onsager, Grystal statistics. I. A two-dimensional model with an order- disorder transition, Phys. 

Rev. (2) 65 (1944), 117-149. 
[38] Y. Peres, Lectures on "Mixing for Markov Chains and Spin Systems" (U. of British Columbia, 

August 2005). Summary available at http://www.stat.berkeley.edu/~peres/ubc.pdf. 
[39] C.-E. Pfister, Large deviations and phase separation in the two-dimensional Ising model, Helv. 

Phys. Acta 64 (1991), no. 7, 953-1054. 
[40] C.-E. Pfister and Y. Velenik, Interface, surface tension and reentrant pinning transition in the 2D 

Ising model. Comm. Math. Phys. 204 (1999), no. 2, 269-312. 
[41] C.-E. Pfister and Y. Velenik, Large deviations and continuum limit in the 2D Ising model, Probab. 

Theory Related Fields 109 (1997), no. 4, 435-506. 
[42] L. Saloff-Coste, Lectures on finite Markov chains. Lectures on probability theory and statistics 

(Saint-Flour, 1996), Lecture Notes in Math., vol. 1665, Springer, Berlin, 1997, pp. 301-413. 
[43] S. B. Shlosman, The droplet in the tube: a case of phase transition in the canonical ensemble. 

Comm. Math. Phys. 125 (1989), no. 1, 81-90. 
[44] H. Spohn, Interface motion in models with stochastic dynamics, J. Statist. Phys. 71 (1993), no. 5-6, 

1081-1132. 

[45] N. Sugimine, A lower bound on the spectral gap of the 3-dimensional stochastic Ising models, J. 

Math. Kyoto Univ. 42 (2002), no. 4, 751-788 (2003). 
[46] L. E. Thomas, Bound on the mass gap for finite volume stochastic Ising models at low temperature. 

Comm. Math. Phys. 126 (1989), no. 1, 1-11. 

E. LUBETZKY 

Microsoft Research, One Microsoft Way, Redmond, WA 98052-6399, USA. 
E-mail address: eyalOmicrosof t . com 

F. Martinelli 

DiPARTiMENTO Di Matematica, Universita Roma Tre, Largo S. Murialdo 1, 00146 Roma, Italia. 
E-mail address: martin@mat.uniroma3.it 

A. Sly 

Microsoft Research, One Microsoft Way, Redmond, WA 98052-6399, USA. 
E-mail address: allansly@microsoft.com 

F.L. Toninelli 

CNRS and ens Lyon, Laboratoire de Physique, 46 Allee dTtalie, 69364 Lyon, France. 
E-mail address: fabio-lucio.toninelli@ens-lyon.fr 



